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Synopsis : 

A new computer code has been written for the NASA Lewis Research Center 
named HEPSPARC (which stands for HEat Pipe SPA ce Radiator Code). This code is 
used for the design and analysis of a radiator that is constructed from a 
pumped fluid loop that transfers heat to the evaporative section of heat 
pipes . 

This manual is designed to familiarize the user with this new code and to 
serve as a reference for its use. This manual documents the work done thus 
far. It is also intended to be the first step towards verification of the 
HEPSPARC code. 

Details are furnished to provide a description of all the requirements 
and variables used in the design and analysis of a combined pumped loop/heat 
pipe radiator system. A description of the subroutines used in the program is 
furnished for those interested in understanding its detailed workings. 
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1.0 GENERAL OVERVIEW: 


1.1 INTRODUCTION: 


The demand for increased spacecraft power brings with it a rise in waste 
heat rejection requirements due to power conversion system inefficiencies. 

This growth of power systems necessitates the ability to analyze the numerous 
radiator possibilities available to meet the requirements. The most efficient 
method to meet the specifications must be determined from the existing 
options. Thus, analytical modeling tools have become a requirement for power 
system heat transfer analysis, optimization, and comparison. 

Existing spacecraft heat transfer systems are pr ima rily composed of 
pumped loop systems (1), heat pipe systems (2), or combinations of the two, 
while advanced, low mass systems are also being developed (3, A, 5) as 
alternatives. Current radiator designs show an increasing dependence on 
systems employing heat pipes. The growing dependency on heat pipes is in part 
due to the nearly isothermal character and high heat transfer capabilities 
inherent in the heat pipe concept. The increased dependence is also due to the 
fact that each heat pipe of a radiator system can be considered an independent 
entity. A designer can plan for heat pipe failures due to environmental 
hazards over the radiator life without concern for overall heat rejection 
system failure. Both these facts allow a heat pipe radiator to typically 
become lighter and smaller than a comparable pumped loop system. 

However, a radiator comprised strictly of heat pipes may have limited 
application due to complex power system/heat pipe interfaces. A combination of 
a pumped loop with heat pipes allows for greater heat rejection system 
versatility. Therefore, a new heat transfer analysis computer code has been 
developed for the NASA Lewis Research Center that analyzes a radiator 
comprised of a pumped fluid within a transport duct that moves across the 
evaporator section of a heat pipe array. Figure 1 depicts such a system. 


1.2 CODE DESCRIPTION/USAGE 

The new code can be used for numerous applications including the 
validation of contractor designs, in-house design and analysis efforts, and 
basic conceptual analyses. It is designed to allow maximum versatility and is 
capable of, but not limited to the analysis of radiator systems for Stirling 
and Brayton engines, thermoelectric elements, and terrestrial, planetary, and 
orbital situations . 

The code’s radiator analysis allows two options. A detailed thermal 
analysis that determines the overall radiator characteristics, or an option 
that determines a minimum mass radiator. In the latter, a specific parameter 
is varied for the optimization. The program then iterates within specified 
limits until the minimum mass configuration is determined. The specific 
working fluid within a heat pipe only allows operation in certain temperature 
ranges and pressures. Therefore, the code has been written in such a way that 
only applicable working fluids can be used at certain temperature levels. When 
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the radiator must produce a large temperature change in the pumped loop fluid 
the code indexes through different heat pipes chosen for the design. The 
temperature of a particular heat pipe is checked against specified limits to 
determine its acceptability at a particular radiator location. If the heat 
pipe is found to be operating at or over a limit, a heat pipe with a different 
working fluid will be chosen. This capability allows the program to be used 
for the design and analysis of radiators that require a large temperature drop 
in the transport duct fluid. A radiator for a Brayton cycle system presents 
such a situation. The optimization option also accounts for multiple heat pipe 
types in its calculations and allows independent optimizing calculations of 
each heat pipe type to minimize the entire radiator mass. 

One of the most important aspects of a design for a space based radiator 
system is the manner in which it is protected from the space environment. The 
protection method used is referred to as armoring or bumpering. The 
determination of the required armoring in the code follows the methodology 
described in Reference (6). It uses an iterative technique to balance the 
number of allowable losses of heat pipes or transport ducts with that 
predicted by an assumed armor thickness using a "critical particle flux". 
Particle impact effects are accounted for by the choice of two models: the 
NASA thick plate, with specified adjustment factors (7,8), or the NASA thin 
plate model (8). Armoring thickness is adjusted until the design specified 
number of heat pipe losses is at its maximum allowable value. This technique 
allows the analyst to account for the radiator environmental conditions. It 
also allows a reduced mass design by use of redundancy in the heat pipes as 
well as with the transport ducts. Additional program options will allow the 
use of offset bumpering as transport duct protection. 


1-3 CODE HIGHLIGHTS /DETAIL; 

The complexity and generality of this program causes its input to be 
detailed. Table 1 summarizes the input and user options available. An 
adequate knowledge of heat pipes is needed to specify all the required 
parameters. Basic understanding of space radiator design is also necessary to 
make maximum use of the program. A significant amount of detail about the 
radiator is returned from the code analysis (see Table 2). When the 
optimization option is used, information about the iterative step values is 
also output in order to track the minimization process. 

The code offers many unique features. A built-in library offers the 
designer a wide selection of radiator materials. The option of using 
properties not already in the library by allowing user defined properties is 
also available. This allows a designer to easily determine the relative impact 
of new materials on a radiator design. Data is also available within the 
program to indicate known compatibility information between the materials 
chosen for use. 

Other program features include the use of fin and heat pipe efficiencies. 
This accounts for the non-isothermal fin conditions as well as radiation 
interchange between the fins and adjacent heat pipes. 
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TABLE 1: INPUTS AND USER OPTIONS 
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SUPPLEMENTAL MASS ADDITIONS 
RADIATOR SPACE CONSTRAINT SITUATIONS 



TABLE 2: LIEPSPARC OUTPUT 



5 


BEGNNNG AND END OF MSSJON CHARACTERISTICS 



Any radiator geometry is easily modeled by using specified input and 
options available. The basic configurations analyzed thus far include a simple 
flat plate design, a cylindrical shape, and a conical shape. Each of these 
configurations still allow for the use of all the program options previously 

mentioned • 

The basic structure of the radiator code's heat pipe analysis section is 
based on portions of the NASA Lewis Heat Pipe Code (9). The technique employed 
for evaluation of the heat pipe characteristics is an iterative method of 
determining the thermal resistance and temperatures of the various sections of 
the heat pipe. Essentially, a one-dimensional, lumped parameter, therma 
network model is used to obtain the heat transfer characteristics of the heat 
pipes. Heat transfer limits are checked using fundamental equations available 

in the literature. 

The Lewis Heat Pipe Code allows for the use of a wide variety of heat 
pipe wick structures. The choices available include wrapped screens, with or 
without arteries, sintered metals, with or without arteries, and grooves. An 
option also exists that allows the user to define all the parameters needed 
for the thermal calculations. This provides the flexibility required to 
evaluate new or different heat pipe concepts. 


1.4 PURPOSE OF THIS MANUAL 

The new code, named HEPSPARC, for HE<U Pipe SPA« Radiator Code, is very 
generic in nature and allows for the design and analysis of a combined heat 
pipe/pumped loop radiator for a wide variety of situations. This user manua 
describes the usage, capabilities, and structure of this code. Details of all 
aspects of the code will be described. Sufficient input requirements and 
specifications will be provided to supply information to allow a new user t 
work with the program. 

This document allows users of various interest levels to pinpoint the 
sections they may wish to read. As already seen. Section 1 has provided basic 
overall description of the code and its application. Basic material involve 
with the workings of the program is furnished in Section 2, while Section 
details the requirements for actual usage of the code. Section 4 provides 
information for understanding the code output as well as additional basic 
information for its usage. Section 5 provides some final instructions and 
guidance for the user. Section 6 contains (in the original version of this 
manual only) a floppy disk with a copy of the program along with a description 
of the subroutines involved in the code. Section 7 contains the is o 
references for this document. 
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2J) CODE DETAILS : 


This section of the user manual describes the basic calculation methods 
employed in HEPSPARC. Detailed equations, methods, and procedures are 
furnished to help describe, document, and justify the workings of this 

at ° r desi 8 n code. Details of the structure, models, and layout of 
HEPSPARC are furnished. 


lii GENERAL: 


As mentioned previously, the program is designed to be generic in nature. 
Therefore, much is left up to the designer and his or her understanding of the 
code. Sufficient flexibility exists to handle a great many situations that 
may arise. However, as with any computer code, the scope of possibilities 
addressed has to be limited. 


2 . 1.1 Operating System 

The code is currently limited to a VAX computer system. This does not 
necessarily mean that it can not be placed on other machines. However, in its 
current form, the code makes use of an IMSL (Integrated Hath Subroutine 
Library) subroutine which is present on the NASA Lewis VAX system. 
Specifically, the function used is entitled DPOID. It evaluates the value of 
a Poisson distribution for the armor thickness calculations. Any transfer of 
the HEPSPARC program to a different system would have to account for this use 
of the IMSL function. Additionally, the use of this program on a personal 
computer is not recommended at this time. Work is still required on the code 
m order to make its operation more efficient. Certain calculations currently 
used in the code make it more adaptive to a mainframe machine. Until work can 
be done to optimize the actual workings of the code, any attempted use on a 
personal computer may prove frustrating to a user. 


2 . 1.2 Basic Procedures 


This code is used for the design of a heat rejection radiator system, 
primarily for space applications comprised of pumped fluid loop that 
transports heat from a heat source to the evaporative section of heat pipes. 
The pumped fluid can be either a liquid or a vapor with the allowance for 
phase change to occur in the loop. The heat pipes are devices that use the 
latent heat of a fluid to transport heat from one point to another. Fluid 
evaporates from the evaporative section of a heat pipe’s wick structure and is 
carried down the pipe by hydrodynamic forces until it reaches the condensing 
section. At that point, the vapor condenses, thereby giving up its latent 
heat which is radiated to space from the outside surfaces of the heat pipes. 
The fluid from the condenser is then returned to the evaporator section via 
capillary action in the wick structure and the process starts over. 
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Heat Transfer 

The heat pipes act as the primary location at which heat is rejected from 
the system. The pumped fluid loop is contained within what is referred to as 
^transport duct or header. This also rejects heat, but its contribution to 
the overall system heat transfer capabilities is typically small when compared 
to that of the heat pipes. For space applications , the heat t^nsfer occurs 
by radiation exchange with the environment. An option g 

convective coefficient for the external surfaces of a radiator. This was 
added to allow the code to model terrestrial applications. 

The code allows for two types of computational design options. The first 
option allows one to specify the radiator design. The second option allows 
for the code to produce a design that has a minimum mass. For the first, all 
the basic characteristics of the radiator are specified. The size of the 
transport duct as well as all the characteristics of the heat pipes are known. 
This type of analysis is used as an analytic tool rather than a design tool. 

It can help determine the performance of an existing radiator design, 
allows itself to be used specifically for the verification of the size 
performance of previously designed radiator systems. 

I„ this case, the heat transfer of the system is specified as well as the 
required temperature change of the transport duct fluid. The of 

this duct fluid establishes the fluids that may be used in the system heat 
pipes. This method of analysis follows the transport duct fluid through the 
radiator. Each heat pipe is analyzed for its heat transfer capabilities The 
amount of heat transferred is then associated with the temperature change 
experienced by the flowing fluid using the simple relation: 

Q - mCpAT 

This process continues with more he.t pipes being edded until '“herthetota 1 
heat transfer of the system is met. or the tmper.ture of the transport duct 
fluid is such that a different heat pipe fluid is required. In the latter 
case the analysis proceeds in the same manner and continues to index throug 
the user specified heat pipe fluids until the required heat transfer of the 
system is met. The calculated number of each type of heat pipe ("type here 
refers to a heat pipe that utilizes a different working fluid) is the minimmn 
required that must survive the entire mission duration. T^ s val ^® 
with the armoring calculations (described below) and specified f® d “" d f^ y 
determine material thicknesses. Variation in the armoring thickness has an 
effect on the heat transfer calculations, and therefore, the process is 
iterative in nature in order to determine the final radiator performance. It 
slolldle noted that the heat transfer of the transport duct itself is also 
included in the determination of the overall system heat transter 
capabilities . 

The second method available is for the design of a radiator to meet 
specific mission needs. With this, one of five possible radiator 
characteristics is to be varied to obtain a minimum mass radiator system 
desien Currently, the variables available for optimization include 
systL’ redundancy (for both the heat pipes as well as the transport ducts). 
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heat pipe diameter (and as this variable changes, the internal thermal- 
hydraulic characteristics of the heat pipe also change), fin width (defined as 
the physical spacing between the heat pipes), and fin thickness. The latter 
three variables can be optimized for each type of heat pipe used in the 
radiator system. To allow for the additional calculations in the optimization 
analysis to occur in a timely manner (i.e. the basic mass comparison of 
several different radiator designs), each individual heat pipe is not 
examined. Rather, a representative heat pipe is chosen for a section of the 
radiator transport duct fluid that undergoes a specified temperature change. 

An estimate of the number of heat pipes needed for the heat transfer of the 
radiator section is determined. As before, this continues until either the 
specified heat transfer of the system is met, or a new heat pipe working fluid 
is required. Also as before, the number of heat pipes determined is used in 
an iterative manner with the armor thickness calculations. The transport duct 
transfer is also included in the optimization analysis. With the 
variables available for optimization, a radiator design can be progressively 
re ^* ne< * or der to obtain the overall optimal design configuration. It may 
be suggested that when this is complete, the simple mass determination mode of 
the program may be used to verify the results from the optimization algorithm. 
It is desired to expand this section of the program substantially. First, the 
addition of heat pipe length as an optimization variable is desired. Then, 
the program could examine all the optimization variables itself to determine 
the final design. This would significantly reduce the amount of user 
interaction required by the program but would prove to be a formidable 
programming task. 


Basic Armoring Methodology 


Radiators in orbit around the Earth are exposed to an environment of 
micro-meteoroids and space debris particles. These objects, though typically 
very small, may be travelling at velocities in excess of 20 km/s. Upon 
contact with a radiator surface, there is often sufficient energy transferred 
to puncture the material. If the puncture occurs at a location that is a 
fluid boundary, the fluid is lost to space and the radiator area that was 
serviced by the fluid becomes ineffective. Therefore, the design of space 
radiators must accommodate the effects of this environment. 

The material used in a radiator system protects it from the space 
environment. When the radiator material is used to protect the underlying 
objects, such as a fluid boundary, it is referred to as armoring or bumpering, 
depending on the application of protection methods (this area will be 
addressed in more detail later in this document). To prevent all punctures of 
a system, the armor must be made very thick. This results in a heavy radiator 
design. Redundancy, as described above, can reduce the mass of a radiator. 
Some components may become inoperative due to the punctures that may occur, 
but the radiator still is able to function due to the excess number of 
components. The result of using redundancy is that the armoring of the 
components becomes significantly less, and an overall radiator system is less 
massive. 


Heat pipes represent components that are generally small in size that 


9 



operate independently from one another. Heat pipes can operate in parallel 
with one another providing redundancy in a system. This redundancy can be 
used both for a pumped loop and heat pipe radiator design. The code allows 
for the design of a radiator that is comprised of a number of independently 
pumped loops, or transport ducts, that are coupled to their own set of heat 
pipes (see Figure 2). Thus, it is possible for the designer to allow the 
transport ducts to have "extras" as well as to have redundant heat pipes 
attached to them. Both of these variables are options that exist for the mass 
minimization design in the code. 


2.1.3 


GENERAL FLOW CHART 


The code is broken down into two main sections. These are the mass 
minimizing use of the code (referred to as the optimizing method) and the 
basic code use (referred to as the non- optimizing method) . These make use of 
a number of common subsections for their calculation purposes. These include 
input and output sections, heat pipe and transport duct heat transfer 
calculations, materials data, armor determination, material compatibility 
checking, stress analysis, and pressure drop calculations. 

The flow chart of the operation of HEPSPARC showing the control structure 
and basic interaction with the various subsections is shown in Figures 3 and 4 
for the optimized and non-optimized cases respectively. More detailed 
information pertaining to calculational methods will follow as each of the 
various program sections is described. The sections are each made up of 
several subroutines. The purpose of each subroutine (listed by name) is given 
in section 6.1 for reference. 


2.1.4 


DESCRIPTION OF CODE WORKINGS 


The following sections describe the logic and course followed in the 
operation of the code. The information given is an expansion of the flow 
chart presented previously. It furnishes information about the intricate 
workings of the code for those who are interested or require this level of 

detail . 


2. 1.4.1 


NON-OPTIMIZED METHOD 


The first step taken in the design/analysis of a combined heat 
pipe/pumped loop radiator is to estimate of the number of heat pipes required 
to transfer the specified amount of heat. This value * 8 based on “V 
life mission requirement that needs to be met and would be part of the basic 
design specification. The total heat transfer required is divided by bh ® 
number of independent transport ducts that are specified to survive for the 
mission. The analysis proceeds, concentrating on the design of one transpor 
duct and its heat pipes and then uses symmetry to determine the overall 
radiator design. 

A representative radiator temperature is chosen and the heat that one 
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FIGURE 2s BASIC DESCRIPTION OF RADIATOR 
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FIGURE 3: NON-OPTIMIZED FLOW CHART 
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FIGURE 4: OPTIMIZED FLOW CHART 
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heat pipe can dissipate is determined. This calculation is done assuming the 
heat pipe has no protective armor. The value calculated is used to estimate 
the total number of heat pipes required for the entire radiator section. If 
the temperature and temperature change required in the radiator design 
requires more than one type of heat pipe working fluid, then the switch-off 
point (location where a different heat pipe fluid is used) is determined. 

This allows the determination of the required heat transfer of each group of 
heat pipes in a radiator section. The same procedure as outlined above is 
used for determining an initial guess for the required number of each heat 
pipe type to meet the specifications in the radiator section. 

Once an initial guess for the number of all the heat pipe types is 
obtained, redundancy is included, and a total area is determined for the 
radiator section. This area is used along with the requirements on mission 
life, system altitude, redundancy, and required probability of radiator 
survival to determine the armor thickness for the heat pipes. At this point, 
a suitable estimate exists about the radiator section to begin the detailed 
analysis. A new value for the number of heat pipes is determined accounting 
for the armor thickness. At this time, an estimate of the required size of 
the transport duct is determined based on the number of heat pipes required. 
With this information, an estimate is made of the required armor thickness on 
the transport duct. This is based on the user specified transport duct design 
(i.e. location of the heat pipes in relation to the transport duct), duct 
material specifications, redundancy, and required survivability. When the 
armor thickness is determined, it is used to calculate a thermal resistance. 
With this known, an estimate of the transport duct heat transfer is calculated 
based on the average duct fluid temperature in the radiator section. This 
value is used to adjust the number of heat pipes required in the system. The 
iteration process with the heat pipes starts again and continues until the 
final design is converged upon for total number of heat pipes and transport 
duct size. 


2. 1.4. 2 OPTIMIZED APPROACH 

The basic approach used for the analysis with the optimizing approach is 
not significantly different than that described above for the basic analysis. 
Several differences do exist though. 

To be able to do the calculations for the radiator design numerous times 
for the mass minimization, groups of heat pipes are used instead of single 
ones to determine heat transfer capabilities. The total temperature drop of 
the radiator pumped loop fluid is broken into increments of approximately 5 K. 
Allowances are made when switch-off points between different heat working 
fluids are used. The required heat transfer of the fluid within this 
increment (i.e. product of mass flow, specific heat, and temperature change of 
the fluid) is determined. The arithmetic average temperature of the fluid 
within one of these increments is used as the temperature at the outside of 
the heat pipe evaporator. The heat transfer capabilities of a single heat 
pipe at that average temperature is calculated. The total number of heat 
pipes for the increment is determined by multiplying the ratio of total 
required heat transfer in the increment to that obtained for a single heat 
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pipe. This process continues for the entire radiator temperature drop and all 
heat pipe fluids. The mass of this radiator design is then calculated in the 
same way as is done for the non-optimizing approach. 

The optimized design approach requires the user to choose upper and lower 
limits for the variable that is changed during the mass minimizing process 
(e.g. a minimum and maximum redundancy value must be specified when this is 
chosen as the optimizing variable). The system mass is initially determined 
for the specified upper and lower values as well as for the mid-point 
optimizing variable value. Brent’s optimizing method determines the final 
variable value to give a minimized system design. This method uses the slope 
of system mass lines to determine the next variable guess for minimizing the 
system mass (see Reference 17 for additional information on this numerical 
method). 


The technique used for minimizing the radiator system mass only allows 
one variable to be varied at a time. Therefore, when more than one heat pipe 
working fluid is required, the optimizing variable of only one heat pipe type 
is used to determine a minimum mass system. When this calculation is 
completed, the variable relating to the next heat pipe type is optimized to 
calculate the new system mass. This continues until all the different heat 
pipes types have been optimized for the specified variable and the final 
result reflects the overall system minimum mass. 

It is not desirable to have the optimizing algorithm terminated due to a 
self imposed error in the program (refer to section 4.3), therefore, certain 
flags and conditions are set internal to the program during the analysis. Any 
situstion which produces an error condition during the normal analysis will 
tri^8®r a flag in the optimization run that establishes a high, artificial 
value of the system mass for that design point. The variable causing the 
error condition is eliminated from the analysis. This continues until a 
optimizing variable range is established that is with the program geometrical 
constraints. Therefore, even though the user may specify the optimized 
radiator to be established within a wide range of variable values, the program 
may reduce this range to one that can be addressed by the limitations of the 
code. In the output, the user is able to follow the events of the program to 
tell what has happened during calculations with out-of-bounds conditions. 

This method has worked a majority of the time, however, several instances have 
been found in which the optimization has not progressed correctly using this 
procedure. Therefore, it is recommended that when an error condition is 
triggered (indicated by an estimated mass of the system of 1E8 kilograms in 
the output), the user should rerun the analysis with a smaller specified range 
for the optimization variable. 


2.2 GEOMETRY POSSIBILITIES 

The radiator configuration can be varied substantially. The only basic 
constraint is that the design is one in which a fluid is transported to the 
evaporative sections of heat pipes. Various geometries can be designed via 
the use of the two view factors prescribed for the radiator system. Anything 
from flat plates, to cylinders, to cones, to a spoked wheel design, as well as 
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others, can be modeled with the use of independent view factors for two sides 
of the radiator. A user must spend some preliminary design effort to 
determine the correct values for the two view factors. 


Within any geometry, the actual radiator operation can be accomplished in 
many ways. As mentioned previously, redundancy can be used for both the heat 
pipes and the transport ducts. Therefore, any number of independent transport 
ducts can be specified in a design, with any number of these being required to 
survive the entire mission. Attached to these ducts are the heat pipes that 
can also be redundant. These redundant heat pipes will simply increase the 
basic size of the transport. 


The value of the heat pipe redundancy is expressed as a percentage of the 
number of heat pipes required for successful system operation at the end of 
mission. For instance, if it is determined that 100 heat pipes are required 
to reject the heat in a radiator section, and 82 redundancy is specified, 
there will actually be 108 heat pipes present in the final design. Due to the 
random nature involved in a failure from particle impacts, redundancy is used 
for all the heat pipes with different working fluids. Therefore, the same 
relative amount of extra heat pipes will be present for all the different 
working fluids. 


2 . 2.1 


TRANSPORT DUCT 


The transport duct carries the fluid heated by the heat source to the 
evaporative section of the heat pipes. The analysis done in the program does 
not assume any specific configuration of the transport duct or its interface 
with the heat source. The user information about the transport duct treats it 
as an arbitrarily shaped channel vith a specified hydraulic diameter and 
wetted perimeter. Figure 5 shows the possible locations of the heat pipes 
relative to the transport duct. The heat pipes may simply be attached to the 
top of the duct or may actually penetrate the duct wall. In the former case, 
the heat pipes are conductively coupled to the transport duct fluid while in 
the latter situation the pipes are coupled by forced convection to the duct 
fluid. Additional permutations exist for the transport duct /heat pipe 
interface. The heat pipes may be attached to the transport duct on one or two 
sides. Refer to Figure 6 for an example of these two cases. 

The final possibility dealing with the transport duct geometry is the 
allowance of a return duct in the system (refer to Figure 7). This situation 
may allow a better representation of an actual radiator geometry require n a 
system. The return duct may or may not have the same shape and size as the 
primary transport duct. Its specifications are entered independently into the 

code . 


2 . 2.2 


HEAT PIPES 


The heat pipe calculations and the geometries of both wick and envelope 
are extracted from the Lewis Heat Pipe Code. This allows for a multitude of 
possibilities in the heat pipe design. The choices of wick design include the 
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FIGURE 5: HEAT PIPE /TRANSPORT DUCT 

CONFIGURATIONS 



(A) HEAT pipes protrude into duct 



FLUD FLOW 


(B) HEAT PIPES ARE LOCATED "ON" DUCT 
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FIGURE 6: HEAT PIPES ON ONE 
AND TWO SIDES OF TRANSPORT DUCT 



(A) TWO SIDES 



(B) ONE SIDE 
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FIGURE 7: TRANSPORT DUCT DESIGN WITH RETURN DUCT 
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RETURN DUCT 



use of axial grooves, wire mesh screens, with or without arteries, sintered 
metal wicks, also with or without arteries, and a generic wick in which the 
user specifies all the geometric and hydrodynamic aspects of the wick 
structure. The basic assumed shape analyzed for the heat pipe envelope is a 
cylinder, however, specifying certain variables can result in the modeling of 
non-circular envelopes for the heat pipe cross section. The details of the 
requirements to describe non-circular heat pipes are described in Section 3.2. 
Also specified by the user are the heat pipe entrainment number, wetting angle 
of the fluid /envelope combination, and the nucleation radius available on the 
heat pipe envelope surface. These factors are also described in more detail 
in the user input sections of this manual. 


2.3 MATERIALS 

A built-in materials library is featured in the code. It i6 comprised of 
a substantial number of materials that have been considered for use in 
radiator systems. The choices available are described below with any 
applicable comments effecting their use. 


2.3.1 FLUIDS 

2. 3. 1.1 TRANSPORT DUCT 

The transport duct fluids available for use in the code are shown in 
Table 3. The code number for the fluid is used to identify it in the program 
input and output. 


Table .3 : Transport Duct Fluids 


Fluid 


Code Number 


Lithium 

NaK 

Helium 

Toluene 

Sodium 

Potassium 

Cesium 

Water 

Mercury 

Methanol 

Freon-11 

Heptane 

Dowtherm A 

Ammonia 

User specified properties 


1 

2 

3 

A 

5 

6 

7 

8 
9 

10 

11 

12 

13 

1A 

15 


When the user input option is chosen, the program requires thermodynamic 
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information of the fluid that is to be used in the calculations The 
- ---- * eeded is ^^ific heat, density, thermal conductivity and 
Jherefire' * emperature de P endenb properties are not available for’use 

ra *S e - d uVpro^rZs^ 

keit^nind J he . analyais may not be as accurate as desirSd! This must^e 
working f luidsZhich i! S1 ? n ‘ J he Program is easily expanded to include more 
current l^t. ot the JrograT ' th ° Se liSted * b °'" the 


2 . 3 . 1.2 


HEAT PIPE 


fi A ™ llable beat Pipe duct fluids are shown in Table 4 . Included in the 
odified to add additional fluids as information becomes available. 


Table 4_ Heat Pipe Fluids 


Code Number 


Potassium 2 

Lithium 2 

Sodium 3 

Cesium 4 

Water 5 

Mercury 5 

Methanol 7 

Freon -11 3 

Heptane 9 

Dowtherm A 20 

Ammonia 1 1 


Applicable Temperature Range 


773 

K 

to 

1273 K 

1273 

K 

to 

2073 K 

873 

K 

to 

1473 K 

723 

K 

to 

1173 K 

303 

K 

to 

473 K 

523 

K 

to 

923 K 

283 

K 

to 

403 K 

233 

K 

to 

393 K 

273 

K 

to 

423 K 

423 

K 

to 

668 K 

213 

K 

to 

373 K 


2,3.2 

2 . 3 . 2.1 


METALS 


LINERS /ARMORS / FINS /WICKS 


The metals currently available for use in the code for liners (the use of 
liners for heat pipes and transport ducts is discussed in Section 23 2 4 ) 

Zbe r r in th^ terial ’/ nd f - S are ShOWn in Table 5 * The — ocia^ed refereZ 
sho^ All o? th t0 ldentdfy them in the Program output and input is also 
K f 86 Mterials be used *°r the transport duct, heat pipes 
including the wicks, or fins. Any of them may be used as a liner or armoring 
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material in the system design ^ever it should ^ noted that the^etals^^ 
with an asterisk do not have known armor property capaoiiici 
coefficient will be determined by the following equations 


Kl - 0.92/(p/1000 .)°' 5 with p - density in kg/m 3 


K, is the impact coefficient for the NASA Thin Model - ref. er to 
section 2 . 5.2 of this manual for detailed discussion of the 
mathematics of the calculations for impacts. 

coefficients for materials with known 

Thi. equation is a curve fit^of ^impact ^coefficient,^^ ^ ^ ductile 

mater iale^however, the u.e o£ the material. »ith uhlcncwn amoring propertie. 
should be carefully scrutinised when used in a radiator design. 

Table 5 j_ Materials for Liners, Armor,. fins , and Wic ks 


Metal 


Code Number 


Stainless steel 
Titanium 
Beryllium 
Aluminum 

Magnesium lithium 
Niobium-1 zirconium 
User properties 
Iron** 

Nickel** 

Copper** 

Tungsten** 

Inconel** 

Molybdenum** 

Tantalum** 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


** 


the program, and currently, the use of these materials do not 
provide any armoring protection- 


The user specified properties option allows the incorporation of a ^y 
material for use in the program. This, therefore, allows the r P 
such things as composite materials for use in a radiator d ®£ign- 

SSsSi" 
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currently available. Therefore, the same warning exists for materials that 

have very temperature dependent properties used in a system that experiences a 
large temperature change. a 


2-3. 2. 2 BUMPERS 


As of yet, no calculations are included in the code to allow for use of 
offset bumpers in the determination of armoring requirements. Input exists in 
the program for the use of bumpers, but the information is not used. Only 
basic armoring calculations are assumed to occur. 


2 . 3 . 2. 3 


COMPATIBILITY INFORMATION 


Much experimental work has been done with heat pipes, and it has been 
determined that not all heat pipe working fluids and heat pipe envelope 
materials are compatible with one another. Two materials may produce an 
interaction that causes the heat pipe to fail. A listing of known 
compatibilities (acceptable and adverse) is contained in Reference 11. This 
information has been included in the HEPSPARC code. Based on the material 
choices, information is provided in the code output to indicate known 
compatibility information between the heat pipe fluid and its wick and 
envelope material. Even though the transport duct and transport fluid do not 
behave in the same manner as the heat pipe, compatibility information is also 
provided for the duct system to alert the designer of any known history 
regarding the materials chosen for the overall design. 


2_. 3 . 2 . 4 LINERS 


A material is chosen for a 
maximum armoring capabilities, 
with the fluid it must contain, 
the system (refer to Figure 8). 
fluid and outside material. It 
with the best heat pipe fluids 
lmer material will also offer 


heat pipe or transport duct system to provide 
However, often this material is not compatible 
In such a case, a liner may be employed in 
This provides a material barrier between the 
allows the use of the best armoring materials 
even though they may be incompatible. The 
some armoring capabilities. 


. k _ T ^ 1S code allows for the design of transport duct system and heat pipes 
with the liner material. Even though the user may specify a particular 
combination of armor, liner, and fluid, it still may not be manufacturable. 
Nothing is contained in the code to estimate the compatibility between the 
armor and liner material. Several things are recommended to the designer to 
evaluate this aspect. One is to examine the coefficients of thermal expansion 
of the two materials and the stress that would be induced in them for the 
radiator design. If one of the materials will yield, the use of the material 
combination may be questionable. An additional area to check is the 
possibility of a chemical interaction between the liner and armor material. 
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2. 3.2.5 


PROPERTIES 


The material properties used for the mlv.i, ^ fki . 

ZTlT d*.:. 00 "' I? 1 '* ° l the " de <«*« to .ection 

(Reference 9). Additional data was the resul^of mater H Plp6 C ° de 
obtained from various literature sources (Referenc^) ial property curve fits 

obtained from a NASA report on the subject (Referetct ei ? WaS 


2.4 

2.4.1 


heat transfer details 

OVERVIEW 


determine thriett 0 Lttsf“ n ^at d occurs e from S** 11 * he “ ethods em P lo y ed to 
to provide sufficient infoLuon and reference: "“tn" T^l 11 
to folio, the program calculations and";:::::::. interested user 

progr^r^r sssKiSs d r ctibea - ibe 

vers d° n ° £ . the code - T »““:Lrr«ut. t ::e D t 

sssi' ~ - sr *“-* 


2.4.2 


2 . 4 . 2.1 


transport duct fluid 
basic equations 


trans^t^tl^^o'::^::^.:: the temperature drop from the 

as-; ~ -ft™: r he 

4 ns usea in the program are described below. 

Cylinder; SS 3 ^ PlP# evaporator: Heat transfer from fluid to a 


Liquid metals: 

h - (k/D outside of heat pipe )[5.o + 0.025(Re D Pr) °- 8 ] 
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Other fluids: 


h - (k/D„ tlla . o, t„t ^HO.me^Pr'") 

.here Re p - Reynolds number based on heat pipe evaporator 
diameter 

Pr - Prandtl number 

Heat transfer from fluid to transport duct: 

Liquid metals: 

- ,wn _ a )f 6.133 + 0.015(Re D Pr)°- 8 ] 

h * (k/D outsid* of hast pip* 'l 


Other fluids: 

Laminar flow: 

h - 3.66(.klD tmZ ,u 9 of h*»t pip*) 

Turbulent flow: 

S - Ck/tws. e... plp .)(0.023R.„“-*Pr»->) 

here, Re„ is based on the transport duct hydraulic 
diameter 

when a condensing fluid is stipulated a constant 
coefficient of 15,000 W/m K l8 “ Se ' ffllts for tu0 .phase fluids in a zero 
involved with f “^""'I'ue U ^typical of condensing coefficients present 

0 ”^^"”' ™n be elsily modified when an acceptable correl.tion is found 
for zero gravity situations. 


24.2.2 AREA MULTIPLIER 

A multiplier may be used for Jf^JT^esi stance* between* the transport duct 
option may be exercised when the thermal variable is used as a 

??uid and 7 a heat pipe in the flui^stream^s ^J^^transfer area. 

simple multiplicative fact order to reduce temperature drops 

It increases the area of heat transfer in “cder to rea P that 

that may occur. The use of this *“°" s ,„' i „'“hat increase heat 

have convective fins attached to t e evap without knowing the 

transfer. This variable allows the use of ^ is ippU c.ble for a 

details of the fin geometry. ** 'srtit may have a low 

Brayton system in which the working J t0 re duce the temperature 

A ” aloe o£ i -° 
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for this variable will allow an analysi 
envelope determining the heat transfer 


s to proceed with the 
area . 


basic heat 


2 . 4.3 


HEAT PIPE 


In situations in which more than one heat- *1. u i , 

system, the switch-over point between fluids ** iV® required for a 

upper bounds of the lower ooeratinc tom governed by the recommended 

a radiator to be designed fir an end ™ ^ Pipe fluid * This allows 
specified final system requirements are °° ndition in ^hich 

exists when all redundant heat cim hL! h SUCh a conf juration 

just able to meet the system heat^ransf ^ destr °y ed » and the radiator is 

as high a temperature as oossihlp ? P - f , uid to a lower one to occur at 

be met without much possibility of thf k! requirements will be able to 

their thermal or hydraulic S t f t £% heat £*!>«• operating beyond one of 
mission configuration, the ^"‘flSld”" thl ' a ” d 0f 

be at a lower temperature and will rwulr. a > switch-over points will 

further into its recommended intended design range P ** WlU ° Perate 


2 . 4 . 3.1 


EVAPORATOR /CONDENSFB 


m xjsssrs —r aa is — 

equations will not be furnished here Tt I.*! d therefore, the specific 
The^heat f ° r “tinnl ol'chf hea^pipe. 

r 1 ^T^Ls-hit 1 ~ra*s: sss. 

the outside of the cond“se “ IT “*" “ 

i^ug ^“uisiie^ :„ h d r lysis n fthe heat p^p*"'* tr".^." 

resistance, q^tp^Ket Sf “ 

SETS’ th" d h“l pl;e iSt “rthe rOU8 S tL VaPOr itaal " d ^ “ ^ condenser 
reversed from that in the th condenser the resistance network is just 

through the condenser Sck SETS' TTyT”’!! th \ eapor/liquid interface, 
the condenser. The resistances foIThTr *5 through the liner and armor in 
separate nodes to accSS fol lhe noTlhn??'"', 1 T tion iS d ° ne at tw0 
through the heat pipe at near th.rLl-hyd™Iic iL^LTESI"" dr ° P 


2 . 4.4 


FINS 


from tie ouIs'id'rsuIlIc^To^he^nl^" ° CCUrS C ° ” di *“°" interchange 
assumed to haye ^pIs PiP Ihelrara 
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FIGURE 9: HEAT PIPE THERMAL RESISTANCE NETWORK 
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used in radiator design to help reduce the mass of the system. The fins 
decrease in temperature along their length due to the heat loss from their 
surface. Typically, the effectiveness of the fin is expressed as a fin 
efficiency. This is the ratio of amount of heat a fin actually transfers to 
the amount of heat that could be transferred from a same sized isothermal 
surface at the temperature of the fin base. Due to the fin/heat pipe 
configuration, not all of the heat radiated from the surface will be directed 
towards the space sink (see Figure 10). Some of the energy that emanates from 
the heat pipe will be intercepted by the fin or adjacent heat pipes, and 
thereby reducing the effectiveness of the heat pipe surface. The fin will 
also radiate energy that is intercepted by the heat pipes. This phenomenon is 
accounted for in the code. An analysis was done to determine the effects of 
geometry on the fin and surface efficiencies. The description and results of 
that analysis are presented below. 


2. 4. 4.1 GRAY BODY ANALYSIS 

The procedure for determination of surface and fin efficiencies is 
outlined in Reference 13 and was performed on the basic configuration shown in 
Figure 10. The surfaces were assumed to be gray, i.e. they were not assumed 
to be perfect emitters and absorbers. The analysis broke the configuration 
down into differential sized elements. View factors were determined between 
all elements, and an iterative scheme was developed to determine the 
radiosities of the differential elements. From that, the total heat 
transferred from the surfaces was determined. This analysis was done for 
various combinations of fin length, thickness, pipe size, and surface 
emissiyities. Additionally, the analysis allowed adjacent heat pipes to be 
operating at slightly different temperatures which may actually be the case in 
a radiator. The heat transfer values calculated were then compared to a gray, 
isothermal surface of the same size to establish the effectiveness of 
the heat pipe, via a surface efficiency, and the fin, via a fin efficiency. 
These values are used in the code in a manner described in the following 
section. 


2-4. 4.2 TABULATED DATA 

The results of the gray body analysis has been entered into HEPSPARC by 
means of a simple look-up table. Geometrical factors, temperatures, thermal 
conductivity, and surface emissivity values are all used to determine the 
appropriate table location in which to obtain surface and fin efficiencies. A 
factor to account for the effect of non-isothermal adjacent heat pipes is also 
in tabulated form, and is used as a multiplicative adjustment to the 
efficiency values . 

The total heat transfer that occurs from the outside of the heat pipe 
surface is calculated from the data obtained in the table as: 

Q “ surface ~ T sink )( *1 surf ace^heat pipe + ^fin^fin) 

with Fj or 2 * view factor of side 1 or 2 of the radiator to 
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FIGURE 10: HEAT PIPE/FIN BASIC ARRANGEMENT 
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SINK 



L sink 


o * Stefan-Boltzman constant 


^surfscs * ^sink 
*)*urf*c« ♦ *1fin 


^hsat pip* » Af in 


heat pipe and sink temperatures 

adjusted surface and fin efficiencies 
respectively. These already account for 
non-isothermal adjacent heat pipes and a 
specified emissivity to be present on the 
surfaces. That is why e does not appear 
in the above equation. 

Total heat pipe and fin surface areas 
respectively 


design uses different fin and base materials. U Wh * n radlator 


2.5 


ARMOR DETAILS 


radiator ° £ *""° r thlck "- s •* the 

described in th^ IS 7 he pipes * n<J ““»P°rt ducts, is 


2 . 5.1 


FLOW CHART 


in Figure 11^ The'detaH th ? calculation of the armor thickness is shown 


2 . 5.2 


£- ASI £ ARMORING equations and methods 


directly ^ 1 e f «e™Y.Vi;!'^ MtlOT ° f th ' *™° r is taken 

this document. Several modijii.%- P 7 t) ' n 0 P sl " d k“e for completeness of 

Micro-meteoroids : 

Particle distribution models exist that describe the level of certain 


31 



FIGURE 11: DETAILED ARMORING CALCULATION 
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mass particles present in the space 
near earth flux has been modeled as 


environment . 


For micro-meteoroids, the 


lo 8ioN t - -14.37 -1.21og 10 m 


with N t * 


average meteoroid flux (number of particles 
to or greater than m grams, per m 2 -s) 


of mass equal 




G * ' °' 568 + °-‘ 32,R '* * <«78 + H,/637 8 „ . orbital altitude in to, 

E “ (1 + cos6)/2 with 0 ■ arcsin[6378/(6378 + H)] 


impacting particle * “ 


K i">p 0 - 35 Pp 0 - 17 V 0 - 8 e 


where 

t « material thickness, cm 

Kj- penetration constant of given material 
nip- mass of particle, g 
P p - density of particle, g/cm 3 
V - particle velocity, km/s 


the^ss^f a A par^iJle P iequired e Jo JJ*® eq “ tlon . c * n be solved to determine 
is substituted into th. ,? d to P e "«tcate a given material thickness. This 
having ^ "" uaTto o StS’SST.i: ««.rmine the flu, of particles (JS”,* 

material thicLess f or Z " "n ^UyT) ' •' ’ **• 

yields: velocity (Vi). Combining all the above 


N 


ti 


4.27 x 10- 15 G.Ef 3 - 45 K, 3 - 45 o °- 5 8v 0.31 

* i Pp v i 


corresponds to a particular velocity increment 
t ° t f 1 J? lue of the fll “. the particle velocity profi 
* ach velocity is then used t0 ~ 


In order to obtain 
is broken apart and 
an overall lethal 


N lt - I N ti Pi where Pi . probability of velocity Vi occurring for 

the micro-meteoroids 


When this is done for 
expression results: 


the micro-meteoroid distribution. 


the following 


N lt - 7.15 x 10‘ n G e Et' 3 - 45 K 1 3 - 45 p p 0 - 58 

The particle density, Pp , is assumed to be 0.5 g/cm 3 . 
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~ ss-ss.". as.i.-as ?ri:i.r.s^.su" 


or fewer lethal 


(h<G) 


- exp(-N lt AT) 


r-n 

I (N lt AT) r /rl 
r-0 


where 

^(h<n) 


Nit 

A 

T 


probability of impact by n or fewer lethal particles ‘ 
corresponds to the probability of survival of the radiator 
heat pipes or transport duct systems, 
lethal particle flux, particles /m -s 
exposed radiator area to be protected, 
mission duration, seconds 


nr 


An iterative technique is used J^^r.isf^'caUut.tion. 

equation. The required *" a ^ d “'™tsumed design (for the transport ducts), 
(in the heat pipe case), or increase this area. The number of 

The redundancy of the system is above equations. For a 

redundant components represents t e va u ^ t ^ e probability of n 

given thickness, redundancy, at u . probability is then compared to 


Debris : 

Space debris is addressed in the same way as the micro-meteoroids, 
total lethal flux for debris is determined to be: 


The 


N lt - 8.89 x lO-Ut- 2 ’ 22 *! 2 ’ 22 ?*' 31 


;; £: 

present for year 2020. 

I£ the radiator exposed areas were the s.»e for ^^"o-neteorcid^.nd 

fSr.urr sse 

described in Section 2*5.5. 

space dehris and oicro-^t.oroid 

P-ic. a conservative 
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design. 


2 . 5.3 


BASIC BUMPERING METHODOLOGY 


structure being protected ! he A U certain r di r t that iS . not direc tly adjacent to the 
underlying surface. The concept behind ® tance . exi ® ts between the bumper and 

bumper will cause impacting pa^icle lo eTtL'rL ^ the U8e ° f the 
with the bumper or to be broken Mr! h be vaporized due to impact 

lower mass and be disperseJ over thlMdl ? " Sulti “8 fra g*«nts will be of a 
thinner than if it were exposed^ t l h Tl 3erl ^ n f * urface - This surface may be 

the particles reaching it Lve a lower .nd^ere^™ ° f 8inCe 

particle flux* mass, and therefore, a lower lethal 

Reference 8, however /as^enti^ned 1 ^^*^ 60 ^!^ 88 ° f bumper is 8 iven in 
part of the code calculations Thic c !• th 18 currentl y unavailable as 
code aware of the need and desire to includ" ?£? added t0 “ ke the user of the 
improved version of the Radiator code ^ arm ° rin8 ‘ Spect in an 


2 . 5.4 


USE OF MULTI-LAYER MATERIALS 


surface l! “» d 

specified between lethal particU flur .^ ^^ relationship 

Therefore, an alternative r a ir,!f, d thickness may be conservative. 

the equations determined above were for the thi!* 8 been established * Whereas 
penetration (NASA Thin Plate Model) ^he^a multi^ 8 ° f “ ate f ial at incipient 
required thickness -7 be Educed ^ lieferenceMK' 1 * 3 " 1 "' ater1 * 1 

sesaJ^inu; y ^o/^A‘r?hic k ^ h Lt:^^ 0 f/ n .t imP,Ct ° f * p * rticle * 

is (Reference 6): P1 * te M ° del) * The basic equation relating these 

P. - K„m 0 - 35 p m 0 - 17 V 0 - 67 

where P„ - depth of penetration, cm 
K. - material constant 
m “ mass of micro-meteoroid, g/cm 3 

?' - impac ^velocity denSit5 '' * V,1 “ e ° f °- 5 * s •»*>«. 

determine C the t wa 11 ^thickness ^when's ???*? ° f P enetr *ti°n is then applied to 

walled structure. niCk " eSS Mhen a niUti-layer material is „ se d for a thin 
t - CFrP. 

micro-meteoroid'distribution? b It°ia dMe^S"^, 1 ' 11 ’* 1 ““ f ° r th6 


35 



„ lt . 5.79 X 10->7 G.Et-’-"CF>-‘ 5 K. s -‘V“ 


For this equation, the particle density. P p . is again assumed to be 
0.5 g/cm 3 . 

the only choice currently available. 

The same procedure used above may be foiled for the space debris flu* to 
give: 


Mu - 2.69 x IQ' 10 t' 2 ' 22 CF 2 ' 22 K. 2 ' 22 Pd°' 37 


Again, the grand total of these tv^fluxes ^ 1 e q* a "^ 6 i “ e ed%of lethal 
flux'determination'i^the SrltSn" calculation, when a correction factor is 
specified . 

The material thicknesses determined 

Thin and Thick Plate models) f any liner on the section being 

thickness to account for the p converted into an equivalent armor 

■nalvzed . The thickness (t) of a liner is convex 
thickness by use of the following relationship: 


t.quival.nt .nnor ” K 1 .amortlinar /Kl.linar 


ss s a - 

This equivalent thickness is subtracted fr^ the calculated required thichness 
value to determine the actual armor thickness. 


2.5.5 


VULNERABLE AREAS 


_ The area used in the determination of e 
calculations may be varied somewha . therefore, a lethal particle may 

considered to be isotropic in nature, end , d V ot the heat 

arrive from any direction. Ther '*°' ' . d fo/determining the armoring 

pipes and transport duct is required to be used for t he exposed 

The area to be used for space J** 1 LTmor/planar^anner sugges^y 
whether debris behaves ls °“°P^ a y the total area would have to be 

some (Reference 15). In the former * w<mld const itute the exposed 

used, but in the latter, a projected a (multiplied by the exposed 

surface . In that case, the space debris flux, 

length) is used for the determination of armor tor v 
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a • A11 mu 6 f reas used are corrected for the configuration of the radiator 
design. The heat transfer view factor of the two "sides" of the radiator 
surface is used to reduce the exposed area. This is done to compensate for 
the fact that a radiator design may help shield itself. For instance, a 
cylindncally shaped radiator will obviously not have the same probability of 
impact on its interior surface as it does on its outside. Therefore, it is 
assumed that the same view factor exists to the debris and micro-meteoroid 

th ^ 6X r StS f ° r the surfaces t0 the specified heat sink temperature. 
Additionally, for the transport duct, some additional shielding will occur if 
a return duct is added to the system. This shielding will reduce the exposed 
area of the duct. The amount of the shielding will be determined by the 
specific configuration of the duct design. Since it would be very difficult 
t0 determine all the data required to calculate the amount of 
shielding that occurs when a return duct is used in a system, it has been 
assumed that 75 Z of the total surface area of a round or non-round transport 
will be exposed to the space environment fluxes of debris and micro-meteoroids 
when a return duct is used and the total surface area is used for armor 
determination. When a projected area is assumed for use with the micro- 
meteoroids, only half of this value of vulnerable area (37.5Z) is used for 
armor determination on the transport duct. These values were determined to be 
present on specific designs HEPSPARC was used with. It is believed to be a 
good estimate of the vulnerable area present when a return duct is used. If 
the user does not feel that it correctly models a particular design, it may 
easily be changed within the program. 


2.6 OTHER DETAILS 

This section details other features of the code. This includes pressure 
drops and other details of the transport duct fluid as well as details 
concerning stress levels calculated in the code. 


2.6.1 PRESSURE DROP 


The pressure drop experienced by the fluid in the transport duct is 
estimated in the program. This allows for the determination of pumping power 
(see below) so a designer may estimate this aspect of the radiator system. 

The pressure drop and pumping power will detract from the overall system 
performance and will have to be supplied by some system power source. 

The pressure drop in the system is dependent on the frictional energy 
loss the fluid undergoes while moving through the transport duct. The code 
breaks the calculation down into several different sections for analysis. 

When the heat pipes do not project into the transport duct, the total pressure 
drop is determined by the basic duct flow equation: 

AP “ 4pV 2 fL/(2D) where V “ fluid velocity 

f ■ friction factor (see below) 

L « length of duct 

D - hydraulic diameter of duct 

p = density of fluid in transport duct 


37 



The equations used for friction factor vary in the above analysis with respect 
to the Reynolds number of the flow. They are determined as follows: 

For Reynolds number (Re) < 5000 

f - 16/Re 

and for 5000 < Re < 30,000 
f - 0.079Re' 0,25 
and for Re > 30,000 

f - 0 . 046Re"°" 20 

When the heat pipes project into the duct, the Reynolds number used is 
d.-nSaM upon the Sri. being examined. areas required to be earned include 
the entrance and exit of the radiator panels, the area in between heat pipes, 
and the area directly adjacent to the heat pipes. A Reynolds number 
determined for each of these sections and the corresponding friction factor 
determined. Each of the pressure drops associated with these areas is 
determined and added to estimate the overall pressure drop exp * ct ®j J* 
entire system. The pressure drop for the total number of panels determines 
the value for each independent radiator loop* 

The friction characteristics for the pressure drop determination were taken 
from Reference 16. Although the correlations are specifically for round 
pipes the use of the hydraulic diameter should give reasonable estimates o 
the friction factor for the pressure drop calculations for non-round transport 

ducts . 

The pressure drop due to heat pipes being in the duct is accounted for by 
determining the pressure drop over a single tube in cross-flow. It is 
realized that this is not the best way to model the situation due to P os ® 2ble 
effects from the proximity of the transport duct. Some type of single, in- 
line tube bank correlation may be better, but none was found for incorporation 
into the code. It is estimated that the use of the single tube in cross-flow 
im produce conaerv.tive results. Tube banks typically have a lower pressure 
drop due to a shielding effect from upstream tubes. 

Loss coefficients can be added to the system pressure drop to account for 
the presence of elbows, restrictions, and other actual duct possibilities. 

The way in which they are used is to multiply the loss coefficients and 
velocity head of the transport duct fluid to obtain the additional system 
pressure drop. No use of the equivalent diameter method is employed. 


2.6.2 PUMP WORK 

Pump work is estimated based on the calculated system pressure drop. It 
is the volume flow rate of fluid (Q) multiplied by the pressure drop: 
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Pump work - QAP « mAP/p with m » 


mass flow rate of fluid in transport 
duct 

P ■ density of fluid in transport duct 
AP* system pressure drop 


^o„!wr teS ^ he POW6r required for * 100Z efficient pump to move the fluid 

any additional' n^nirJ l **' designed ’ Actual requirements will depend on 
y additional pumping requirements present in the loop as well as the 

electric * 1 efficiency of the pump end motor combinttion choeen 


6 • 3 STRESS ANALYSIS 


The hoop stresses in heat pipes and transport ducts are estimated in »-h» 
program. Even though the transport duct, and at times the heat pipes may not 

the *ff Ul ? r Shape ’ the stress determined in a cylinder is used to estimate 
the effect of pressure in the system. estimate 


Hoop stress: a - Pr/t with P - pressure of fluid in transport duct 

or heat pipe 

r * radius of heat pipe (based on 

diameter of heat pipe used for armor 
determination) or hydraulic diameter 
of transport duct divided by 2 

* thickness of wall of transport duct 
or heat pipe (includes liner and 

armor thickness) 
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3.0 USER INPUT 


3.1 


INTRODUCTION 


This section is used to provide a line by line description of the how the 
code is used. Both the interactive and non- interactive input characteristics 

will be explained. 

The non-interactive method will be addressed initially. The procedure 
that will be followed will be to show the actual lines of data ret ^j- r «* 
the nroeram The meaning of each variable will be explained as well as a y 
aoolicable variables or flags that may be set or triggered within the program 
b? P a given response. This will help in the description of the J n £ era ^* JJ®* 
of the program. All combinations of responses will be addressed for both the 
optimization and non-optimized use of th. coda. The ™ •• they appa.r 

from the code will appear in bold faced text as will the lines of data that 

are created in the input file. 


32 INTERACTIVE AND NON- INTERACTIVE METHODS 


3.2.1 


NON- INTERACTIVE USE 


This section describes the actual physical make-up of the input file. 

This file can be created manually, but the interactive use of the pr°g«* 
llso accomplish the task in an error free manner. The details in this section 
wiU ill™ the user to easily modify a file created during interactive usage 
of the urogram or actually create a file from scratch. It is recommended 
that the interactive method be used at least for the initial ^t^section 
radiator design to ensure the data has been input correctly. 

an^ input f inline (shown in bold face) will be expressed with the *PP™P^** e 
variable names, then the options, values, and ramification of certain choices 
will be described as they relate to the variables. This will serve two 
purposes . First, it will compliment the information provided below concerning 
the "walk through" of the program. Second, it will furnish a list of 
«ri.bUs th« are .cruelly used in th. coda for c.lcul.tional *'>'“081' 
purposes! This -ill aid a user in the det.min.tion of the logic : involved in 
Lmv of the algorithms. Two additional things need to be specified. The code 

uses an unformatted read statement of the input file, ther ® f °^®’ * 
not have to worry about ensuring that certain variables are located in 
specific columns of the file. The input information must be separated y 
rorrani or a snace. Secondly, the code uses an implicit specification for 
integers and^eal valued variables. Therefore, any variable beginning wit 
the letters I through N (inclusive) are considered to be integers and must be 
entered into the input file as such. These variables generally constitute 
counters and flags within the program structure. 


HOPT , NPAN , PWID , IPANT , DIH , HPERIM , ICOND , M , COEF 
NOPT - 


Optimization flag variable. NOPT - 0 for non-optimized use of code 
while NOPT - 1 is to use the algorithm in a mass optimization 
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manner . 


NPAN - 
PWID - 


I PANT - 
DIH - 


Number of radiator panels per independent transport duct path. 

Maximum tota i radiator panel width in meters. Can only be used 
when NOPT-O. Occurs when a fixed transport duct width will be 
specified for the radiator (IDES-1; explained below) and this 
value will be used to determine the heat pipe condenser length. 
This value will depend on if the heat pipes are specified to 
emanate from one or two sides of the transport duct. 

Type of radiator panel. IPANT-1: heat pipes project from only one 

- ?! transport duct; IPANT-2: heat pipes project from both 
sides of the transport duct. 


Primary transport duct effective hydraulic diameter in meters. It 
is determined from the following equation: 


HPERIM - 
ICOND - 


M - 
COEF - 


DIH - 4A C /P W where A c - duct cross-sectional area 

P w m duct wetted perimeter 

Wetted perimeter of primary transport duct in meters. 

Fluid phase flag. If ICOND - 0: transport duct fluid is single 
phase (i.e. all liquid or all gas); ICOND - 1: two-phase fluid is 
used in the transport duct. Currently, when this option is 
exercised, a heat transfer coefficient of 15,000 W/m 2 K is used for 
the condensing vapor. This is due to the fact that a zero-gravity 
correlation condensing correlation is unknown to the author. 

Number of different heat pipe types used in the radiator. A heat 
pipe type is defined as one in which the working fluid is 
different. The code is intended to allow the use of several (up 
to 5) different heat pipe working fluids within on radiator system 
m order to allow a large transport duct fluid temperature change. 

Loss coefficient of an independent transport duct network (non- 
dimensional) for supplemental pressure drop determination. 


If the optimization routine is used (N0PT-1), the following lines follow in 
the input data set: 


NVAR 


This variable indicates which variable is to be used for the mass 
optimization. The available choices are: 

NVAR - 1; Fin thickness 

NVAR « 2; Fin width (defined as the distance between outside 
diameter of heat pipes) 

NVAR * 3; Heat pipe redundancy 

NVAR * 4; Heat pipe diameter 
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NVAR - 5; Variation of number of independent transport ducts 
surviving the entire mission duration 

If NVAR - 1, the next lines in the input file will be: 

TFINHI , TFINLO 

These are the upper (TFINHI) and lower (TFINLO) constraints of the fin 
thickness (in meters) for the optimization. The algorithm w on y 
determine a minimum mass system that has the optimization variable 
between^hese values . The line is repeated H (as defined above) times 
to account for all the different heat pipe types. The variable order is 
from the highest temperature heat pipe working fluid to t e owes 

If NVAR - 2, then the next lines will appear: 

FWIDHI « FWIDLO 

These are the upper (FWIDHI) and lower (FWIDLO) constraints for the fin 
width (in meters). As with the fin thickness, the line is repeated M 
times starting with the high temperature heat pipe conditions. 

If NVAR ■ 3, then the following will appear: 

REDHI.REDLO 

These are the upper (REDHI) and lower (REDLO) constraints of the heat 
p^reLdancy" It is assumed that all the different •‘e.t ptpe types 
will have the same relative amount of redundancy. The values that ar 
to be entered are calculated as follows: 


RED * 1.0 + Redundancy (in percent) / 100. 


Therefore, a 5 percent redundancy constraint would be specified by a 
1.05 value entered into the input file. 


iiuic t a. * 

value entered into the input file. 
If NVAR ■ 4, the following is specified: 


DHI.DLO 


These are the upper (DHI ) and lower (DLO) constraints of the actual heat 
pip^diameter ^meters. As with the other per. meters, it is tepe.ted M 
times account for the different heat pipe types. It starts with the 
high temperature pipe and continues to the low. 


If NVAR - 5, no additional constraints need be specified. The number of 

=«ihi^^ 

appear in the output file. 

After the lines for the optimization case (if they are required), the 
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following will appears 

HL.ENLT , EXLT , IHLM , HLT , ICONV , LHPLOC 


HL - 


ENLT - 
EXLT - 

IHLM - 

HLT - 
ICONV - 


Maximum radiator duct length per panel, in meters. To be used 
when NOPT-O. Occurs when a fixed transport duct length will be 
specified for the radiator (IDES-1; explained below) and this 
length will be used to determine the heat pipe spacing. The inlet 
and outlet lengths of the transport duct are subtracted from this 
length for each panel, and the heat pipes are then evenly 
distributed over the remaining duct length. If the IDES-1 option 
is not exercised, the required transport duct length will be 
determined for the system based on specified fin width values. 


Entrance length of transport duct for each radiator panel that 
contains no heat pipes (in meters). 

Exit length of transport duct for each radiator panel that 
contains no heat pipes (in meters). The entry and exit lengths of 
the transport ducts may represent portions of a radiator that are 
required for some other use such as a deployment mechanism. 

Transport duct liner material code flag. See 2. 3. 2.1 for the 
associated material flag values. If no liner is to be used, any 
integer value may be entered* 

Transport duct liner material thickness in meters. When no liner 
is used, a value of 0.0 should be specified for the variable. 

Flag that specifies the boundary condition on the external 
radiator surface. ICONV - 0: represents a condition of pure 
radiative heat transfer from the outside surface; ICONV-1: a 
combined convective and radiative heat transfer system exists on 
the radiator surfaces. 


PLOC - Flag that specifies the relative location of the heat pipes on the 
transport duct. LHPLOC-1: the heat pipes protrude into the 
transport duct and come into physical contact with the duct fluid; 
LHPLOC-2: the heat pipe is located on the transport duct and is 
coupled to the duct fluid by conduction through the duct and heat 
pipe wall. This latter arrangement is the one present on the SP- 
100 radiator configuration. 

If ICONV - 1, the next line will appear: 


HCONV 


This is the external convective heat transfer coefficient that is 
present on the radiator external surfaces. 


IHAR, IHWF, AQ , TIN , TOUT , IDES , NPARA , NLOSS , NPAN 
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IHAR - 
IHWF - 
AQ - 


Transport duct material code flag. See section 2. 3. 2.1 for the 
actual material associated with each variable value. 

Transport duct working fluid code flag. Section 2. 3. 1.1 describes 
the fluids and their codes available for use in the program. 

Desired heat transfer of the radiator at the end of life condition 
in kilowatts. 


TIN - 
TOUT - 
IDES - 


NPARA - 
NLOSS - 

NPAN - 


Inlet temperature of the transport duct fluid in K. 

Desired outlet temperature of the transport duct fluid in K. 

Flag that specifies how a radiator design will be developed. This 
option is not available when the optimization method of analysis 
is chosen (i.e. NOPT - 2) and the IDES flag should then be set to 
be equal to 2. When the non-optimized version of the code is 
used, this variable may be set to equal to 1. When this is t e 
case it is assumed that severe constraints are to be placed upon 
the maximum size of the radiator structure. Therefore, both a 
maximum transport duct length (per panel) and maximum radiator 
panel width is assumed to be required in the design. These 
stipulations then allow the program to determine certain variable 
values that would normally be specified as input by the user when 
this option is not used. Specifically, the heat pipe condenser 
length is calculated as the specified panel width (PWID) minus the 
transport duct diameter. This length would be divided by two if 
the heat pipes project from both sides of the transport duct. 

Also, the fin width is determined by subtracting the outside 
diameters of all the system heat pipes and the lengths of the 
inlet and outlet sections of the transport ducts from the total 
transport duct length (accounting for all the radiator panels) and 
dividing this value by the total number of heat pipes "qu^ed in 
the system. Thus, it is assumed the heat pipes and fin will 
up all the available duct area. If IDES is set to equal 2 .both 
the fin width and heat pipe condenser length must be specified by 
the user for each heat pipe type. When IDES-1, all the different 
heat pipe types will be assumed to have the same condenser length. 
This may create some out of limit conditions for the heat pipes. 
The use of the IDES-1 option does have some applications in a 
design situation, but it is recommended that the program be used 
with the IDES variable set to equal 2. 


Number of independent parallel transport duct paths in the 
radiator at the beginning of the mission. 


lumber of allowable inoperative (due to micro-meteoroid and space 
lebris impact punctures) transport ducts at the end-of-mission 
:onf ieuration. 


Number of sections or panels of radiator per each independent 
parallel path. 
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0 


If the transport duct fluid was specified to have 

i.e. an IHWF of 15, the following line will next 
file: 


user defined 
appear in the 


properties , 
input data 


CPHF , CONDHF , VISCHF , RHOHL 

CPHF - Specific heat of the transport duct fluid in J/kgK. 

Thermal conductivity of the transport duct fluid in W/mK. 
Viscosity of the transport duct fluid in N-s/m 2 . 

Density of the transport duct fluid in kg/m 3 . 


CONDHF - 
VISCHF - 
RHOHL - 


The following line will then occur either after the NPAN or RHOHL value: 
RPR.TSINK, FI , F2 , Das , FLOW, ITRANS 


RPR - 
TSINK - 


Inlet pressure of the 
used to determine the 


transport duct fluid in the radiator. It is 
pressure stresses of the transport duct. 


Effective sink temperature the radiator panels experience. To 
yield a conservatively sized radiator design, this variable should 
be the maximum value that will be experienced during a mission or 


FI - 


F2 - 


EMIS - 


FLOW - 


Effective view factor of side number 1 of the radiator to the 
effective sink temperature. For the radiator being designed, if 
the heat pipes are located on (versus in) a transport duct, the 
wo si es of the radiator will not be symmetrical. If thi6 is the 
case, si e 1 is defined as the side of the radiator panel having 
the transport duct located on it. It should be noted, when the 

heat pipes are located in the transport duct, this distinction is 
not applicable. 


View factor of side 2 of the radiator to the effective sink 
temperature. This, along with the FI variable allow for the 
design of a radiator in nearly any geometry. 

The emissivity of the external surfaces of the radiator. All the 

material that comprises the radiator is assumed to have the same 
emissivity. 


This variable represents the mass flow rate of fluid through the 
entire transport duct system (i.e. all independent paths) at the 
beginning of the radiator mission life. This only applies when 

rTrnS Uld it n t “ n !? ort duct is specified as being condensed, 

tie pfnu • if ! f Uid in the trans P ort duct is single phase, 
the FLOW variable should be entered as zero. The flow rate is 

specified in kg/s of fluid. The heat transfer of the system is 
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then determined by: 

Q - mhfg 


jji m mass flow 

h fg - enthalpy of vaporization 

«■“* • *^ur srsrs »***--. 

it will have to be protected with armoring or bumpering. The main 
.rl assumed to be l«.t.d .di.c»t to on, another 
in the system design. Therefore, they will tend to shield each 
other from micro-meteoroids and space debris. The act' 1 ® 1 amount 
of shielding depends on the shapes of the return and main 
transport duct. A value of 25 Z shielding has been assumed for 
ducts P of circular and non-circular cross sections. Therefore, 75Z 
of the total combined area of the normal and return duct are used 
for the vulnerable area calculation. This is for the assumption 
of total vulnerable area available for space debris exposure. 

When a projected area is assumed for this calculation, on . y 
Uf this value of vulnerable area is used a rnjor determination 

on the transport duct for the micro-meteoroid fldx. 

If a return duct is specified for the system (ITRANS - 1). then the following 

should be entered on the next data line: 

RETDIH , RETPER 

RETDIH - Return transport duct hydraulic diameter in meters. 

RETPER - Return transport duct wetted perimeter in meters. 

The following line will eppe.t efter either the det. for the return tr.neport 

duct or after the line ending with the ITRANS variable. 


SURVEP , SURVTD .RED , ALT , IENV , DURR , ADWPP , ADSW 


SURVHP 


SURVTD 


This specifies the required survivability of the hut pipes in the 
system^ This value is used in the determination of th« . armor 
thickness of the heat pipe. A typical value would be *?? *° 
represent a 952 probability of the heat pipes survivability. 

This specifies the required survivability of the transport duct(s) 
In the radiator system. This value is used in the determination 

of the armor thickness of the ducts. The vardab J e “ p^e 

the same way as the SURVHP variable. The product of the h PP 
and transport duct survivabilities will yield an overall system 
survivability value. 
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RED 


ALT - 
IENV 

DURR 

ADWPP 

ADSW • 

ADWPS, 

ADWPS 

NARH - 

NAREA 


This variable specifies the redundancy to be used for the heat 
pipes. This variable is also used in the determination of the 
heat pipe armor thickness. The value entered is determined in the 
same way the REDHI and REDLO are determined, i.e. a RED value of 
1.05 signifies a 52 desired redundancy. It should be noted that 
for this variable, as well as any other variable that can be 
optimized, may have any value in the main data deck when it has 
been chosen for use in the optimization. When the optimizing 
option is used, the controlling variable are the limits specified 

for the optimization specified in the third line of the input 
deck. r 


This is the altitude in kilometers at which the radiator is to 
operate during its mission. This value helps determine the micro- 
meteoroid flux that will be experienced during a mission. 


^ 1 TirMf, riable specifies the environment for the radiator system. 

It IEN V - 1, then only micro-meteoroids are assumed to be present 
and the radiator armoring thickness will be calculated based only* 
on this flux. If IENV « 2, space debris will also be considered 
to be present in the radiator environment* 

This represents the desired mission life of the system in years. 

T v 1S ™y iable is Used t0 add wei 8 ht to the system. It represents 
the additional mass (in kg) to be added in the system analysis for 
each radiator panel present in the entire system. This may be 
indicative of some structural mass requirement for each panel such 
as a deployment mechanism required for each panel. 


This variable is also used to account for additional masses 
present in the system that may not be covered by this analysis. 
This variable represents the extra, lump sum, mass (in kilograms) 
to be added for the entire radiator system. 


NARN , NAREA , AMPMTD , NMOD , AMULT 

This variable is also used to add mass to the system. It 
represents the additional mass (in kilograms) to be added during 
the analysis for each independent section of the radiator. It may 
represent, for example, the mass of a pump that would be required 
for each independent transport duct. 

This variable represents the choice of armoring design for the 
transport duct. If NARM - 1, then normal armoring methods are 
used, while NARM - 2 specifies that offset bumpering techniques be 
employed. Currently, no bumpering calculations exist, and 
therefore, NARM should always be specified to be equal to 1. 

This flag represents the vulnerable area determination to be used 
m the armor determination for space debris protection. A NAREA 
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value of 1 signifies the actual circumferential area is to be used 
while a NAREA of 2 indicates the projected surface area will be 
used in determining the vulnerable area of the heat piP« or 
transport duct. This option was added due to the uncertainty 
involved in the nature of space debris. 

ampmtd - This variable also adds mass to the system analysis. It 

AMPMT represents the addition of mass (in kilograms) for each meter of 

transport duct length for each transport duct. It may represent 
some type of mass adjustment, such as a bumper, for the entire, 
unknown transport duct length. 

NMOD - This variable is used to specify the model to be used for the 

armor determination. If NMOD - 1. the NASA Thin model will, be 
used while NMOD - 2 indicates the NASA Thick model with 
correction factor will be used. Section 2.5.2 of this manual 
discusses the details of these two models. 

AMULT - If the NASA Thick model is used, the AMULT variable *P*^f 

correction factor to be used with it. Currently, the NASA Thick 
model calculations have a built-in correction value J- 5 * 
Therefore, the specification of AMULT does not affect it. This 
will be changed in the future. The recommended range of this 
value will be 1.5 to 2.5. If the NASA Thin plate model is 
specified (NMOD-1), then this variable should be set to be zero. 


If the material chosen for the transport duct liner is to be user specified 
(IHLM -7), the next line will follow. 

CONHLM , RHOHLM , HLK1 

Transport duct liner material thermal conductivity (W/mK) . 
Transport duct liner material density (kg/m 3 ). 

Transport duct liner material impact coefficient (as would be used 
with the NASA Thin plate model). 


CONHLM - 
RHOHLM - 
HLK1 - 


If the material chosen for the transport duct armor is to be user specified 
(IHAM *7), the next line will follow. 

CONHAM , RHOHAM , HAK1 

CONHAM - Transport duct armor material thermal conductivity (W/mK) . 
Transport duct armor material density (kg/m 3 ) 


RHOHAM 
HAK1 - 


Transport duct armor material impact coefficient (as vould be used 
with the NASA Thin plate model). 
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? th * **• to be 

d.t. starting with the highest temperstire ^"klng fl^ h.« %£* 

KVAP , DHPA . HPLTE . HPLTC , IHPF , IHPL , IHPA , TOT , PINT , FACTOR 


EVAP - 


DHPA - 


HPLTE - 


HPLTC - 


IHPF - 


IHPL - 


IHPA - 


IWT - 


FINT 


of 1 the’hest 1 pije . th This"is h thi n iensth*of 0 fh the ( * va >’° r * tor Motion 
the transport d'uct depending 'p^e deslgl." ^ in “ °" 

^assuming thl It" 

nO r iiner is to he need, then Ihis IIluT^uS “ 

I° r Il"'r 18 t0 be “ Sed ' the " ^is value should be set equ.l'to 

^.’'"lilao'n'rm-n"!' 10 ^ U " d in 'he particular heat 
y y ecmon 4.3.I.2 of this manual indimt-oc f^ a ri • j 

fi«t*hiV°* USe End th6ir a PP licable temperature ranges 8 The 

£sm5sl-- 

used, this should be set equal to be zero. 

Section*! * p ® c J- f the arm °r material choice for the heat pipe 

2. 3. 2.1 indicates the choices available for the materials. 

fon S rt Variab i e - identlfieS the wick desi 8 n f ot the heat pine 
following choices are available: P P 


The 


IWT 

IWT 

IWT 

IWT 

IWT 

IWT 


1 

2 

1 

1 

1 

1 


Rectangular grooves 
Screen wick 
Sintered metal wick 

Generic wick described by user properties 
Screen wick with arteries 
Sintered metal wick with arteries 


pipe If h :H^ Ckl ^ S i^ n meters ) °f the fin attached to the heat 
pipe. If the fin thickness is to be optimized, then this vlrilus 
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FACTOR - 


should be set equal to zero. 

This variable represents a multiplicative factor to be used on the 
IhU . hpat transfer area on the heat pipe evaporator section. 

5 evaporator ££ SSTSTuS 

2JSTT2 

SSU! Then, fin material may be desired to be added to the 

“”J“ic”*^rir.^nMtiS e i. t d.n«d! a value of on. should 
be specified for FACTOR. 


WM.DHP.IFIN 
TWM - 


DHP 


This variable specifies the wick material flag. Sect f°* 2 ;**™ 
indicatesthe material associated with the various values of XHM. 

This variable specifies the diameter (in meters) of the heat 

S Wit 

variable! bl^for non-round designs, they may be substantially 
different. 

ris s arirr-ru"*;'". ... 

data deck. 


Section 


IFIN - 


FIHW.HPCL 
FINW - 

HPCL - 


This variable is the fin width (in meters). Fin width is defined 
to be the distance between adjacent heat pipes (i.e. ou s 
surface to outside surface). 

This variable represents the heat pipe condenser length (in 

mo tore ^ . 


If the heat pipe liner material is to be of a user 
7), then the following line is required. 


specified material, 


(IHPL - 


CONDL , HPRHOL , HPLK1 

CONDL - Heat pipe liner material thermal conductivity (W/mK) . 

HPRHOL - Heat pipe liner material denaity <kg/m ! ). 


HPLK1 - 


Heat pipe liner material impact coefficient (as would be uaed with 
the NASA Thin plate model). 
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If the heat pipe armor material is to be of a user specified material, 
7), then the following line is required. 


(IHPA - 


CONDA , HPRHOA , HPAK1 


CONDA - 
HPRHOA - 
HPAK1 . 


Heat pipe armor material thermal conductivity (W/mK). 

Heat pipe armor material density (kg/m 3 ). 

Heat pipe armor material impact coefficient (as would be used 
the NASA Thin plate model). 


bhe . fin material is t0 be of from user specified properties (IFIN -7) 
following line will have to appear: * 


with 

the 


CONDF.RHOFIN 

CONDF - Heat pipe fin material thermal 

RHOFIN - Heat pipe fin material density 


conductivity (W/mK). 
(kg/m 3 ) . 


The following lines describe the wick details for the particular heat pipe 
design. Additional details about the required inputs can be obtained from the 
authors of Reference 9. This section of the code was taken directly from that 


If IWT « 1 (Rectangular Groove), one line is required to describe the wick: 
DEEP , WIDGRV , NUM 


DEEP 


Depth of the rectangular grooves (in meters). 


WIDGRV - 


Width of the rectangular grooves (in meters). 


Number of rectangular grooves around the periphery of the heat 
pipe. Note, the number and size of grooves can not exceed the 
available circumference length (i.e. NUM x WIDGRV < v x (DHP - 
2 (DEEP ) ) ) . 


If IWT « 2 (Screen Wick), one line is required to describe the wick: 

DIAM, WIDE, CRIMP, NUMWRP 

DIAM - Screen wick wire diameter (in meters). 

Center to center spacing between wick wires (in meters). 
Crimping factor of screen wire wick. 


WIDE - 
CRIMP - 
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NUMWRP 


Number of wraps of screen specified for the wick design. 


If IWT - 3 (Sintered Metal), one line is also required for the data; 
TWICR.PORW.RSS 

TWICK - Thickness of the sintered metal wick (in meters). 

PORW - Porosity of the wick structure. A PORW value of 0.50 indicates a 

502 porous wick structure. 

RSS - Average radius of the sintered particles which comprise the wick 

structure (in meters). 


If IWT - 4 

(User Input Wick Properties) then five lines of data are required: 

ARE, ARC, PE 

ARE - 

This variable represents the product of the evaporator wick 
permeability and the wick area (m*). 

ARC - 

This variable represents the product of the condenser wick 
permeability and the wick area (m*). 

PE - 

This is the value of the inner perimeter of the heat pipe in the 
evaporator section (for thermal /hydraulic calculations not for 
armor determination) in meters. 

PORE.PORC, 

AREE 

PORE - 

This variable is the porosity of the wick structure in the heat 
pipe evaporator section. It should be entered into the data deck 
in decimal form (i.e. 0.50 * 502 porosity). 

PORC - 

This variable is the porosity of the wick structure in the heat 
pipe condenser section. It should be entered into the data deck 
in decimal form (i.e. 0.50 - 502 porosity). 

AREE - 

This represents the evaporator effective wick thermal conductivity 
(W/m 2 R). 


AREC.AVE.AVC 

AKEC - This represents the effective wick thermal conductivity in the 

condenser section of the heat pipe (W/m 2 K). 

_ This is the effective vapor area available in the heat pipe 

evaporator section (in m 2 ). 


52 



AVC - 


This is the effective vapor area available in the heat 
condenser section (in m 2 ). 


pipe 


RCM, DCRITE , RATAR 


RCM - 
DCRITE - 
RATAR - 


This variable is the maximum meniscus radius available in the heat 
pipe wick structure (in meters). 


This is the value of the critical entrainment 
heat pipe evaporator wick (in meters). 


dimension in the 


This variable represents the ratio of the inside area of the 
evaporator to its open area . 


AREAVE , AREAVC 


AREWE - 


This represents the actual cross sectional area of the wick 

nnroc^ re the ^aporator of the heat pipe. It, along with the 
porosity are used to determine the mass of the wick in the heat 
pipe evaporator region. 


AREAWC 


This represents the actual cross sectional area of the wick 
structure in the condenser section of the heat pipe. It, alone 
with th e porosity are used to determine the mass of the wick in 
the heat pipe condenser section. 


When IWT - 5 
first is the 


(Screen wick with arteries), two lines of data are required, 
same required for the screen wick alone: 


The 


DIAM.WIDE , CRIMP , KUMWRP 


The second line describes the arteries: 


RA.AWT.NA 


Radius of arteries in the wick (in meters). 
Thickness of the artery wall (in meters). 
Number of arteries in the wick structure. 


If IWT m 6 (Sintered Metal wick with 
required. The first one is the same 

TWICE, PORW.RSS 


arteries), two lines of data are also 
as the sintered metal wick: 


The second line describes the arteries: 
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RA, WICCON, NA 


RA - 

WICCON - 


HA - 


Radius of arteries in the wick (in meters). 


Thermal conductivity experimental constant. For metal or felt 
fiber, a WICCON value of 0.34 should be specified. If sintere 
particles are used, a value of 0.53 is to be specified. 


Number of arteries in the wick structure. 


The final line required for the heat pipe definition is then 


RH , THETA , VEBNUM 


RN - 


Minimum 

meters) 


nucleation radius of heat pipe envelope material (in 
A suggested value of this variable is 1x10 meters. 


THETA 


This represents the wetting angle between the heat pipe fluid and 
wick material (in degrees). 


WEBNUM - 


This represents the Weber number at the onset of entrainment for 
the particular heat pipe. 


This completes the data required for analyzing or designing a radiator. If 
additional heat pipes are required for the design, the data on the line 
starting with the EVAP variable would next be entered. The data wou 
continue until the heat pipe fluid is capable of operating at the radiator 

outlet temperature. 


3.2.2 


INTERACTIVE USAGE 


This section describes the questions that are used to query the user for 
a response during interactive use of the program. Even though it is 

recommended that this approach be used (at least in ^ ia ^ } U ai*data deck 
file for an analysis, it was decided to first describe the actual data deck 
that is required for the program before the actual interactive d^etails are 
provided (refer to section 3.2.1). This section is meant to com P 1 ““* n ^ 
previous one. The actual questions that appear during the use of Jj h ® V P ^? 
will appear in bold faced type. Then, the variable being triggered by 
response will be identified as will any results of the Particular variable 
choice. The description of the variables discussed can be obtained from 

section 3.2.1. 

Both the optimized and non-optimized version of the code 
usage and interactive input requirements will be addressed in th 
The optimized approach will be addressed initially. 
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IS IBS BATA FOR IBIS ROB CURRENTLY TH AN EXISTING BATA FILE 7 (Y/N) 

Answering -N- to this question .utcMtic.ll, invokes the 
interactive usage of the code. By re.ponding wlth" -5' the 
program then asks for the data file name *v» * the 

that file for the analysis. ’ then P roceeds to use 


OPTIMIZATION CHOICES: (ENTER 1 OR 2) 

1. MINIMUM MASS OPTIMIZATION RUN 

2. UNOPTIMIZED, SIMPLE MASS DETERMINATION RUN 
ENTER 1 or 2: 


Variable triggered - NOPT 


If optimization is chosen the following will appear 


DESIRED SYSTEM PARAMETER TO BE VARIED: 

1. FIN THICKNESS 

2. FIN WIDTH 

3. SYSTEM REDUNDANCY 

4. HEAT PIPE DIAMETER(S) 

5. NUMBER OF PANELS SURVIVING 
ENTER CHOICE: 


Variable triggered - NVAR 


SYSTEM PARAMETERS: 

TOTAL HEAT REJECTION REQUIRED (kW): 


Variable triggered - AQ 

RADIATOR INLET TEMPERATURE (K): 

Variable triggered - TIN 

RADIATOR OUTLET TEMPERATURE (K) : 

Variable triggered - TOUT 

RADIATOR INLET PRESSURE (kPa): 

(FOR STRESS CALCULATION PURPOSES) 

Variable triggered - RPR 

IS CONDENSATION PRESENT IN THE TRANSPORT DUCT (Y/N) ? 
Variable triggered - ICOND 

If condensation is present, the following will appear: 


55 



CONDENSING FLUID AVAILABLE (kg/s): 


TOTAL SYSTEM MASS FLOW RATE OF 
(i.e. BEGINNING OF LIFE) 

Variable triggered - FLOW 


SINK TEMPERATURE (K): 

Variable triggered - TSINK 

SS. , t3?i < ^S 8 SiS 0 i I S SIDE WITH THE TRANSPORT DUCT WHEN THE 
HEAT PIPES ARE ON, NOT IN, THE DUCT 
ENTER VALUE 

Variable triggered - FI 
VIEW FACTOR OF SIDE 2 TO SINK: 

Variable triggered - F2 
EMISSIVITY OF EXTERNAL SURFACES: 

Variable triRRered - EMIS 


SURVIVABILITY INFORMATION: 

HEAT PIPE SURVIVABILITY PROBABILITY (Z): 

Not^-^i^the^ interactive^ percent vine i. entered. 

TRANSPORT DUCT SURVIVABILITY PROBABILITY (Z): 

Variable triggered - SURVTD 

If redundancy is to be optimized, the following will appear. 

LOWER BOUND FOR REDUNDANCY IN ANALYSIS (I): 

Variable triggered - REDLO 
UPPER BOUND FOR REDUNDANCY IN ANALYSIS (I): 

Variable triggered - REDHI 

If redundancy is not to be optimised . the following will appear instead: 
PERCENT REDUNDANCY DESIRED: 


Variable triggered - RED 

Note - this is entered in percent also 

EARTH ORBIT ALTITUDE (km): 
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Variable triggered - ALT 

ENVIRONMENT: 

1. MICRO-METEOROIDS 

2. MICRO-METEOROIDS AND SPACE DEBRIS 
ENTER CHOICE: 


Variable triggered - IENV 
MISSION DURATION (YEARS): 


Variable triggered - DURR 


PENETRATION MODEL TO USE: 

1. NASA THIN MODEL 

2. NASA THICK MODEL WITH AN ADJUSTMENT FACTOR 
ENTER CHOICE: 

Variable triggered - NMOD 

If the NASA Thick model is chosen, the next line will appear: 

DESIRED MULTIPLICATIVE FACTOR (RECOMMENDED RANGE 1.5 TO 2.5): 
Variable triggered - AMULT 

AREA USED TO CALCULATE VULNERABLE AREA FOR SPACE DEBRIS IMPACT* 

1. ACTUAL, CIRCUMFERENTIAL AREA 

2. PROJECTED AREA 
ENTER CHOICE: 


Variable triggered - NAREA 

IS THERE CONVECTION ON THE RADIATOR EXTERNAL SURFACES (Y/N) ? 
Variable triggered - ICONV 

If there is to be convection, the following line will appear: 

ENTER CONVECTIVE HEAT TRANSFER COEFFICIENT (W/m**2 K) 

Variable triggered - HCONV 

TRANSPORT DUCT INFORMATION: 

TYPE OF TRANSPORT DUCT: 

1. HEAT PIPES PROJECT TO ONLY ONE SIDE 
2. HEAT PIPES PROJECT IN TWO DIRECTIONS 
ENTER 1 or 2: 

Variable triggered - IPANT 
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TRANSPORT DUCT /HEAT PIPE CONFIGURATION: 

1. HEAT PIPE PROTRUDES INTO TRANSPORT DUCT 

2. HEAT PIPE IS LOCATED ON TRANSPORT DUCT 
ENTER 1 or 2: 

Variable triggered - LHPLOC 

IS THERE TO BE A RETURN DUCT ON THE TRANSPORT DUCT (Y/N) 7 
Variable triggered - ITRANS 

T RAN SPORT DUCT PROTECTION SCHWE: 

1. NORMAL ARMORING PROCEDURES 

2. OFFSET BUMPERING IS TO BE USED 
ENTER CHOICE: 

Variable triggered - NARM 

Note - only choice #1 is currently allowed 

TOTAL NUMBER OF PARALLEL PATHS: 

Variable triggered - NPARA 

If the number of transport ducts surviving is not to be optimized, the 
following line will appear. 

NUMBER OF ALLOWABLE PATH LOSSES 

Variable triggered - NLOSS 

TRANSPORT DUCT(S) EFFECTIVE HYDRAULIC DIAMETER (m) : 

Variable triggered - DIH 

TRANSPORT DUCT WETTED PERIMETER (m): 

Variable triggered - HPERIM 

If a return duct is to be used in the analysis, the following two lines will 
appear : 

RETURN TRANSPORT DUCT(S) EFFECTIVE HYDRAULIC DIAMETER (m): 

Variable triggered - RETDIH 
RETURN TRANSPORT DUCT WETTED PERIMETER (m) : 

Variable triggered - RETPER 

ESTIMATE OF TRANSPORT DUCT(S) LOSS COEFFICIENTS (PER PANEL): 

Variable triggered - COEF 
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PANEL INFORMATION: 


NUMBER OF PANELS IN RADIATOR (PER PARALLEL TRANSPORT DUCT PATH): 
Variable triggered - NPAN 

versi °” o£ the co< *® i* »«ing used, the following line 


1 . 


2 . 


DESIGN OF RADIATOR TO BE BASED ON COMPUTER GENERATED 
(i.e. FIXED TRANSPORT DUCT LENGTH w/ NO OPTIMIZATION) 
USER SPECIFIED FIN WIDTH * 


TUBE SPACING 


Variable triggered - IDES 

If IDES - 1 is specified, the following two lines will appear 

MAXIMUM PANEL WIDTH (m) : 

(TO BE USED TO DETERMINE CONDENSER LENGTH) 

Variable triggered - PWID 

MAXIMUM TRANSPORT DUCT LENGTH PER PANEL (m) : 

Variable triggered - HL 

INLET LENGTH OF TRANSPORT DUCT(S) (m): 

Variable triggered - ENLT 

OUTLET LENGTH OF TRANSPORT DUCT(S) (m): 

Variable triggered - EXLT 

TRANSPORT DUCT LINER THICKNESS (m) : 


Variable triggered - HLT 

If there is to be a liner, the following will appear 
TRANSPORT DUCT LINER MATERIAL: 


NOTE: 


CHOICES 8 THROUGH 
DEBRIS PROTECTION 


14 WILL OFFER NO STRENGTH FOR MICRO-METEOROID/ 
IN CALCULATIONS 


A list of material choices will then appear, 
entering the number of the material chosen. 


The user should respond by 


Variable triggered - IHLM 


If the material choice is 
lines will appear: 


the user specified material, the following three 
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TH ERMAL CONDUCTIVITY (V/mK) : 

Variable triggered - CONHLM 

DENSITY (kg/m3): 

Variable triggered - RHOHLM 
T HRESH OLD PENETRATION CONSTANT: 

Variable triggered - HLK1 
TRANSPORT DUCT ARMOR /BUMPER MATERIAL: 

A list of materials will then appear. The user should respond by entering the 
number of the material chosen. 

Variable triggered - IHAR 

If the material choice is the user specified material, the following three 
lines will appear: 

THERMAL CONDUCTIVITY (W/mK): 

Variable triggered - CONHAM 

DENSITY (kg/m3): 

Variable triggered - RHOHAM 

TH RE SHOLD PENETRATION CONSTANT: 

Variable triggered - HALK1 

TRANSPORT DUCT WORKING FLUID: 

A list of fluid choices will then appear. 

Variable triggered - IHWF 

If the fluid choice is the user specified material, the following four lines 
will appear: 

DENSITY (kg/m3): 

Variable triggered - RHOHL 

THERMAL CONDUCTIVITY (W/mK): 

Variable triggered - CONDHF 

FLUID VISCOSITY (Ns/m**2): 
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Variable triggered - VISCHF 
FLUID SPECIFIC HEAT (J/kg): 

Variable triggered - CPHF 
HEAT PIPE INFORMATION: 


t t 4 ^ i i S i POint ’ a11 the h6at pipe informat i°n is going to be determined, 
initialiy, a menu will appear indicating the available fluids that may be used 
for the heat pipes in the temperature range of interest. Also appearing will 
be the range of temperature over which the various heat pipe fluids are 

!?£ f? r example ’ if the inlet temperature of the transport duct is 

tjou K., the following menu will appear: 


POTASSIUM (RANGE 773 K TO 1273 K) 
CESIUM (RANGE 723 K TO 1173 K) 
MERCURY (RANGE 523 K TO 923 K) 


The user should specify the fluid choice 
fluid, e.g. enter A for cesium. 


by entering the number of the desired 


Variable triggered - IHPF 

If the heat pipe diameter is to be optimized during the analysis, the 
following two lines will appear: 

LONER BOUND OF ACTUAL HEAT PIPE INSIDE DIAMETER (FOR OPTIMIZATION) (m) : 
Variable triggered - DLO 

UPPER BOUND OF ACTUAL HEAT PIPE INSIDE DIAMETER (FOR OPTIMIZATION) (m) : 
Variable triggered - DHI 

The following note also appears for this optimization. It instructs the user 
to enter all the additional data concerning this heat pipe as it pertains to 
the lower bounds of the optimization. In other words, if it is desired to 
optimize a heat pipe between 0.005 m and 0.05 m diameter, then all the wick 
details and other heat pipe characteristics should reflect the 0.005 m 
diameter design. 

NOTE: THE OPTIMIZATION BASED ON HEAT PIPE DIAMETER USES A SIMPLE 
RATIOING TECHNIQUE. THE FOLLOWING PARAMETER VALUES ARE TO BE THOSE 
ASSOCIATED WITH THE LOWER OPTIMIZATION POINT HEAT PIPE INSIDE DIAMETER 

lilt TT1? * lm 


The heat pipe diameter optimization has a direct effect of the wick design. 

As the code currently stands, the only type of heat pipe wick currently 
available for optimization is the rectangular groove design. Additional work 
needs to be done to add the other heat pipe wick characteristics to the 
optimization routine. 
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If the heat pipe diameter is not to be optimized, the following line will 
appear: 

ACTUAL HEAT PIPE INSIDE DIAMETER (m): 

Variable triggered - DHPA 

The following additional information is then required for each heat pipe: 

HEAT PIPE INSIDE DIAMETER FOR WICK THERMAL /HYDRAULIC CALCULATIONS (m) s 
Variable triggered - DHP 

LINER MATERIAL: _ . 

HEAT PIPE LINER MATERIAL THICKNESS IN EVAPORATOR (m) : 

Variable triggered - HPLTE 

HEAT PIPE LINER THICKNESS IN CONDENSER (■) : 

Variable triggered - HPLTC 

If there is to be a liner, then the following lines will appear: 

HEAT PIPE LINER MATERIAL: 

NOTE: CHOICES 8 THROUGH 14 WILL OFFER NO STRENGTH FOR MICRO-METEOROID/ 
DEBRIS PROTECTION IN CALCULATIONS 

A list of material choices will then appear. The user should respond by 
entering the number of the material chosen. 

Variable triggered - IHPL 

If the material choice is the user specified material, the following three 
lines will appear: 

THERMAL CONDUCTIVITY (W/mK): 

Variable triggered - CONDL 

DENSITY (kg/m3): 

Variable triggered - HPRHOL 

THRESHOLD PENETRATION CONSTANT: 

Variable triggered - HPLK1 

Then, the armoring information will be obtained: 
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ARMOR MATERIAL: 


A list of material choices will then appear, 
entering the number of the material chosen. 


The user should respond by 


Variable triggered - IHPA 


If the armor material choice 
three lines will appear: 


is the user specified material, the 


following 


THERMAL CONDUCTIVITY (W/mK): 

Variable triggered - CONDA 
DENSITY (kg/m3): 

Variable triggered - HPRHOA 
THRESHOLD PENETRATION CONSTANT: 

Variable triggered - HPAK1 

If the fin thickness is to be optimized, the following two lines appear next 
LOWER BOUND OF FIN THICKNESS (FOR OPTIMIZATION) (m) : 

Variable triggered - TFINLO 
UPPER BOUND OF FIN THICKNESS (FOR OPTIMIZATION) (m) : 

Variable triggered - TFINHI 

If fin thickness is not being optimized, the following line will appear: 

FIN INFORMATION: 

WHAT IS THE FIN THICKNESS (m) : 

Variable triggered - FINT 

FIN MATERIAL: 


A list of material choices will then appear 
entering the number of the material chosen. 


The user should respond by 


Variable triggered - IFIN 

If the fin material choice is the user 
lines will appear: 


specified material. 


the following two 


THERMAL CONDUCTIVITY (W/mK): 

Variable triggered - CONDF 
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DENSITY (kg/m3): 

Variable triggered - RHOFIN 

If the fin width is to be optimised, the following two line, will eppe.rt 

PIN WIDTH IS DEFIHED I* THIS CASE TO BS TH* DISTAHCE BETWEEH ADJACENT 
HEAT PIPES. 

LOWER BOUND OF FIN WIDTH (FOR OPTIMIZATION) (■>: 

Variable being triggered - FWIDLO 

UPPER BOUND OF FIN WIDTH (FOR OPTIMIZATION) (m) : 

Variable being triggered - FWIDHI 

If fin width is not being optimised, the following will sppe.r. 

FIN WIDTH IS DEFINED IN THIS CASE TO BE THE DISTANCE BETWEEN ADJACENT 
HEAT PIPES. 

ENTER FIN WIDTH (m) : 

Variable being triggered - FINW 

. i .-vs r«n ho Bnocified (IDES-1), then the following will 

If the condenser length can be specinea. 

appear : 

CONDENSER LENGTH (m) : 

Variable being triggered - HPCL 

The heat pipe wick end evaporator characteristics will now be determined. 
EVAPORATOR LENGTH IN, OR ON TRANSPORT DUCT («)* 

Variable being triggered - EVAP 
EVAPORATOR AREA MULTIPLICATIVE FACTOR: 

Variable being triggered - FACTOR 

WICK TYPES: 

1 . GROOVE 

2. SCREEN WICK 

3 . SINTERED METAL 

4. USER INPUT WICK CHARACTERISTICS 

5 * SCREEN WITH ARTERIES 

6 . SINTERED METAL WITH ARTERIES 
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ENTER CHOICE: 


Variable being triggered - IWT 

The next thing that will appear will depend on the wick design chosen: 

If the rectangular groove design is chosen, the following three lines will 
appear: 

RECTANGULAR GROOVE DEPTH (m) : 

Variable being triggered - DEEP 
RECTANGULAR GROOVE WIDTH (m) : 

Variable being triggered - WIDGRV 
NUMBER OF WICK GROOVES: 

Variable being triggered - NUM 

If either the screen wick or the screen wick with arteries is chosen, the 
following four lines will appear: 

SCREEN WICK WIRE DIAMETER (m) : 

Variable being triggered - DIAM 

SPACE BETWEEN WICK WIRES (m) : 

Variable being triggered - WIDE 

CRIMPING FACTOR: 

Variable being triggered - CRIMP 
NUMBER OF SCREEN WRAPS: 

Variable being triggered - NUMWRP 

If the screen wick with arteries is specified, the following three additional 
lines will also appear: 

ARTERY RADIUS (m) : 

Variable being triggered - RA 

ARTERY WALL THICKNESS (m) : 

Variable being triggered - AWT 
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NUMBER OF ARTERIES: 


Variable being triggered - NA 


If either the sintered metal wick or the sintered metal wick with arteries is 
chosen, the following three lines will appear: 

THICKNESS OF WICK (m) : 

Variable being triggered - TWICK 

POROSITY OF WICK (Z): 

Note, in the interactive usage of the program, this should be entered as 
a percent value, i.e. a value of 50 would indicate a 50Z porous wick 
structure . 


Variable being triggered - PORW 
RADIUS OF WICK PARTICLES (m) : 

Variable being triggered - RSS 

If the sintered metal wick with arteries is specified, the following three 
additional lines will also appear: 

ARTERY RADIUS (m) : 

Variable being triggered - RA 

NUMBER OF ARTERIES: 

Variable being triggered - NA 

THERMAL CONDUCTIVITY EXPERIMENTAL CONSTANT 

METAL FELT OR FIBER C - 0.34 
SINTERED POWDERS C - 0.53 

ENTER VALUE: 

Variable being triggered - WICCON 

If the wick design is to be derived from user specified properties, then the 
following lines will appear: 

EVAPORATOR WICK PERMEABILITY * WICK AREA (m**4): 

Variable being triggered - AKE 

CONDENSER WICK PERMEABILITY * WICK AREA (m**4): 
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Variable being triggered - AKC 
EVAPORATOR INNER PERIMETER (m) : 

Variable being triggered - PE 
EVAPORATOR WICK POROSITY (Z): 


Variable being triggered - PORE 
CONDENSER WICK POROSITY (Z): 

Variable being triggered - PORC 


s ™ n * s * 

EVAPORATOR EFFECTIVE WICK THERMAL CONDUCTIVITY (W/m**2 K): 
Variable being triggered - AKEE 
COTOENSER EFFECTIVE VICE THERMAL COMBOCTIVITr W«**2 E) = 
Variable being triggered - AKEC 
EVAPORATOR EFFECTIVE VAPOR AREA (m**2): 


They will be 


Variable being triggered - AVE 


CONDENSER EFFECTIVE VAPOR AREA (m**2): 
Variable being triggered - AVC 
MAXIMUM MENISCUS RADIUS (m) : 


Variable being triggered - RCM 
EVAPORATOR ENTRAINMENT CRITICAL DIMENSION (m)s 
Variable being triggered - DCRITE 
EVAPORATOR INSIDE AREA / OPEN AREA: 

Variable being triggered - RATAR 
WICK CROSS SECTIONAL AREA IN EVAPORATOR 
Variable being triggered - AREAWE 
WICK CROSS SECTIONAL AREA IN CONDENSER (m**2): 
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Variable being triggered - AREAWC 


Then, some basic 
will appears 


information required for all the 


different wick structures 


MINIMUM NDCLEATION RADIUS (m) : 
DEFAULT VALUE - 1.0E-6 


Variable being triggered - RN 


WETTING ANGLE: 

Variable being triggered - THETA 
WEBER NUMBER AT ONSET OF ENTRAINMENT* 


Variable being triggered - WEBNUM 


SKSrS ths'hs.t Pipe fluid’s^owest applicable temperature is gc*nt,r^ 

than the desired radiator outlet^ temper. £ th , heat pipe 

fluids will appear, and after th , , design, an 850 K inlet 

design and analysis. For example. ^ the above de *8^^ pipg fluid< lf a 

temperature was specified and c *? desired, the cesium will not be able 

radiator outlet temperature o It8 low gst applicable temperature is 

concerning the design 

of the second heat pipe will be obtained. 

If this second fluid’s temperature range is : shown 

radiator outlet t^erature. then the^ystemm.^ ^ ^ new fluid 

below will appear. Otherwise, the P r ? c * chosen to cover the entire 

choices and will continue until ^.'^“tSd. . a maximum of 

radiator temperature range. As the c P £ a rad iator design. The 

five different heat pipe fluids may be spe Therefore, it 

fluids available cover a temperature range of 213 K to ZO/3 o£ 

is anticipated all radiator designs can be accommodated -^th^ ^ 

fluids. The algorithm was set “P th ^® a.* ir ,n v reauire a large temperature 
radiators for Brayton cycles which will typically require xarg 

change of the transport duct fluid. 


SYSTEM MASS INFORMATION: 

ADDITIONAL MASS PER PANEL TO BE ADDED (kg): 


Variable triggered - ADWPP 
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ADDITIONAL MASS PER METER OF TRANSPORT DUCT(S) (kg/m): 
Variable triggered - AMPMTD 


ADDITIONAL MASS PER SECTION (PARALLEL PATH) TO BE ADDED (kg): 
Variable triggered - ADWPS 
ADDITIONAL SYSTBi MASS TO BE ADDED (kg): 

Variable triggered - ADSW 

The following lines are then used to determine in what file the data lust 
entered will be stored. J 


THIS DATA WILL BE WRITTEN TO A FILE 

WHAT DO YOU WISH THIS NEW FILE TO BE r.AT.i im 7 


If for some reason, an error is encountered in the specified data during 
calculations, the user will still have the data entered stored in a file. 
Next, the user is prompted to determine which file shall receive the output 


NAME OF FILE TO PRINT OUTPUT TO: 


At this time, the computational procedure begins. If an error occurs 
recommended that the output file be examined to determine how far the 
progressed. The next section describes the code output. 


it is 
program 
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4.0 CODE OUTPUT 


4.1 INTRODUCTION 

This section of the user’s manual will describe the output that can be 
expected from using the code. All the flags and descriptors are explained in 
detail, ^example is furnished to aid in this task. It is representative of 
what can be obtained from using this analysis tool with both the optimizing 
analysis and non-optimizing version of the code. 


4.2 DESCRIPTION 

Very little difference exists between the output obtained from using the 
code in an optimized and non-optimized configuration. The difference wi 
seen at the very beginning of the output. The non-optimized output begins 
ilth is«rred(***) heading while the optimization begins with the results of 
calculational steps made in order to determine the system minimum mass. This 
can be seen on the attached example. In the optimized version of the code, 
prior to the starred header, lines will exist that show the optimisation 
variable call NVAR. This indicates which variable was used to obtain «J e 
minimum mass radiator design. The variable value is the same as defined above 
in the input deck description. 

After the NVAR variable is the actual value of the variable being 
optimized. The last entry on the line is then the mass determined for the 
system based on the variable value used in the analysis. Additional U "* 8 
similar to this follow which show similar results for different values of the 
optimizing variable. The mass shown should be decreasing as the optimize 
design isMfined. It should be noted, the calculated radiator mass may in 
fact 8 increase slightly as the optimization calculations proceed. This * 
nuance that was discovered using the version of Brent ’• 

Reference 17. It is due to the convergence criteria that is internaliy 
“pecified in the code to control the extent of the optimixetlon Therefore 
the total absolute minimum mass radiator design may not result, but one that 
is within less than 12 of the minimum value will be obtained. If it 18 ®**? 
that the minimum mass system was not obtained, enough 1 ” £or, ““° n ** *” 
during the optimisation analysis to allow the user to choose the correct value 
for analysis with non-optimised version of the code. It should be noted that 
the last printed value of the optimised variable is used as the basis of the 
entire radiator design and all the subsequent information printed. 

After all the optimizing analysis printout is complete, a starred header 
will appear indicating the total results of the analysis just completed. All 
the information printed from this point on, aside from one small section, will 
then be the same for the optimizing and non- optimizing versions of the co . 

H ^he optimizing version is used, the only difference in the output from the 
starred heading on) will appear at the very beginning. Some statements wi 
appear that indicate the optimizing version of the code was used, and it will 
indicate the variable that was optimized as well as the user imposed limits 

the variable chosen. 
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For the remaining information, user specified requirements and input will 
be intermixed with code calculated values in order to give a description of 
the resultant system design. The first (non- optimized use) or second 
(optimized use) thing to appear will be the basic radiator specifications. 
These include mission requirements and radiator design choices. For example, 
the temperature change and heat transfer requirements of the radiator are 
printed as are view factors, survivability information, and the details of the 
radiator design. Next are the details about the transports ducts including 
materials and fluids chosen as well as basic sizes. Additionally, the 
relative position of the heat pipes and transport duct as well as 
compatibility information is furnished here. 

After the transport duct information appears, all the heat pipe 
information is provided. This includes the basic specified design as well as 
the results of the number of the heat pipes required for each panel in the 
radiator. Additionally, masses and thicknesses are furnished. After all this 
has been provided for the heat pipes, the overall system results are 
furnished. This includes the breakdown of the heat transfer and mass 
associated with each section of the radiator as well as the entire system. 
Additionally, since radiators are compared based on specific mass (kg/m^) and 
specific power (kW/kg), these values are also printed. Several different 
specific area values are provided due to the many different ways in which this 
variable can be defined. These are furnished as system and heat pipe values 
in order to provide an estimate of the contribution of the various components 
to the system mass and design. 

The final information printed is a calculated value of equivalent 
temperatures (based on an area, sink temperature, and heat transfer) as well 
as actual surface temperatures experienced by the heat pipes. The estimate of 
the pressure stresses in the system is provided. 

It is felt that the output provides a substantial amount of information 
for the radiator or system designer to use in the development of a concept 
that could be carried on to a more detailed TRAYSIS or SINDA level of 
analysis. 


4.3 ERROR MESSAGES 

It is hoped that no errors will occur when this program is being used. 
However, certain constraints do exist that cause some out of tolerance 
conditions to occur. These will cause the program to terminate its operation. 
The following are error conditions currently built into the code: 

1) Violation of heat pipe limits 


Heat pipes may only operate within certain prescribed limits for the 
fluids and wick structures described. Therefore, circumstances may allow 
conditions to indicate a heat pipe limit has been violated. If the program is 
being used in the non-optimized mode a statement will be printed out 
indicating which heat pipe limit has been violated on which heat pipe type. 

If an optimization is being performed, an artificial mass (value of 1E8) will 
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be specified for the radiator condition and the program will continue to run. 
By only allowing the choice of the heat pipe working fluids within the 
temperature ranges specified in Section 2. 3. 1.2, the possibility of a heat 
pipe violating one of its limits is substantially reduced. 

2) Material incompatibilities 

This flag does not affect the execution of the program. It is simply an 
output line to flag the designer of a possible problem with the materials 
chosen for the analysis. 


3) Problem with fin efficiency parameters 

The gray body analysis done for the determination of the fin and heat 
pipe efficiencies relies on non-dimensional parameters that depend on the fin 
and surface size characteristics as well as their temperature. The data 
entered into the program to furnish this information is limited. Therefore, 
it is possible that the size and temperature of the heat pipe being designed 
will fall outside the tabulated data range. If this occurs during the non- 
optimized use of the program, statements will appear indicating that the Fin 
or Surface Efficiency Data is out of Tabulated Range”. The user will have to 
change some of the fin variables to rerun the program. The range of allowable 
non-dimensional values used is as follows: 

0.75 < 6 < 0.99 

0.0 < L 2 CTT 3 /k fin t fin < 4 
0 . 025< R/L < 20 

where € - material emissivity 

L - distance between adjacent heat pipes’ outer surfaces 
a “ Stefan-Boltzman constant 

T ■ outside temperature of heat pipe condenser section 
k fin « fin material thermal conductivity 
t fin * fin thickness 
R - heat pipe radius 

As with the heat pipe violations, during the optimization version of the code, 
this error will only cause an artificial high mass to be specified in the 
analysis. It has been determined that an optimal fin size should be well 
within these variables for any analysis performed. 


4.4 PROGRAM APPEARANCE DURING EXECUTION 

Certain other features have also been built into the code to allow the 
user to track the calculations as the program runs on the computer. This, 
first, ensures that the code is operating properly, and second, allows user 
termination of the analysis if something is incorrect. The non-optimizing and 
optimizing versions of the code will appear different (on the terminal screen) 
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i " £0r “ ti °" 18 to track th. 

numhpr^nf £«% n ? n “ opt d version of the code is being used, the actual 
nWntL fj P1P6S calculated to ® eet the system requirements will be 

?? *5* screen \ The user should be to follow the progress of the 

avstei" £?* germination of the actual required number of heat pipes for the 
system. With each iteration, the user should see the number of heat pipes 

closing in on a particular value. If this is not seen within 8 to 10 

ahon^ 1 r S ’ ?! P ^° gram f hould be terminated and the input deck examined. It 

1*1 r !i , ^ ^ W °f k With the program haa indicates the possibility of 

the calculations flip-flopping between two values of the program solution that 
are very near each other due to the fact that the number of Seat pipes h« tJ 

aLSri^ anS \ nUmb6r VSlUe ‘ ThiS haS been addre ssed in the program 

w ? , ? U n0t appear t0 the U8er to occur ®°te than two tLes 

ing the calculations. If it does, the program should be terminated. 

nrin t -!5 en p £? gram f S J used in the optimizing mode, more information is 
Jf! heat Thls ^eludes the number of the heat pipe section being analyzed 
and the heat transfer and temperature of the heat pipe section. However what 
the designer should watch for is the value of the Lss calculated^ the’ 
program. It should be decreasing throughout the program calculations If it 
does not or if the artificial value of 1E8 does Sot change, the usSr wiU 

SSlculStioSS 86 the range ° f the optimizin 8 variable for the program 
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5.0 BASIC INFORMATION. POINTERS, AND GUIDANCE 


5.0.1 


BASIC INFORMATION 


This section is intended to provide some basic information and guidance 
for use of this program. 

The code is currently on the NASA Lewis VAX system. It is hoped to have 
a floppy disk copy of the program attached (see section 6.0) to the original 
of this document in order to make the program accessible to more users. * f 
this is not the case. K. Baker of NASA Lewis has access to the program for 

those interested. 

Note as the code currently stands in the author’s personal directory, it 
is named OPT not HEPSPARC for those who wish to simply copy it directly on e 
NASA VAX system. 

In the event a floppy disk copy is available it needs to uploaded to ‘ 
VAX svstem with an IMSL subroutine library available. At that time, or if it 
is desired to modify this program for some specific use, it will have to e 
re-compiled with the following commands: 

FORTRAN /CONTINUATIONS-50 HEPSPARC* FOR 

The continuations specification is required due to the presence of the 
•look-up" tables for the fin and heat pipe surface efficiencies. It must be 
linked with the IMSL information: 

LINK HEPSPARC, IMSL 

And finally, it may then be run by using: 

RUN HEPSPARC 


5.0.2 


POINTERS AND GUIDANCE 


As mentioned previously, it is highly recommended th.t .^'’.‘design 
created with the interactive usage of the program. It is anticipat g 

will progress by a systematic variation of variables associated wit 
baseline^design^ Therefore, if the date deck is initially created . :hen it 
can be easily! progressively, modified to close in on a final design. Also, 
it is suggested to use the optimizing option of the program initia y 
narrow the design possibilities. Then, for the detailed analysis, the non- 
optimizing version of the code can he used to obtain a *lightlymorerefined 
analysis. Experience has shown this method not only saves design time, it 
also saves computer time. This last fact is especially true when very large 
radiator systems (on the order of megawatts) are being designed and analyzed. 

One last point to be made is some advice on how to interpret the v “ 1 °“ s 
parameters associated with a radiator design and how they are * efe *“ d *° . 
this code. This code was initially written to help analyze and validate the 
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2 st*vs: TJiizruTu' 52 tLiz? v **• 

IVZV^lrX i detr “ ent to the use of tllls 'O' 1 ', and this mnual 
and analysis ’of aloupfad'lu^d’wCf tk ' *“*«" 


hi CONCLUDING REMARKS 

1 £*££ ZttSZ SB 

the 11 .SoulTta Cs^r* ' ’ , ° rkin88 ° f COd ' •" d in Particular ho. 

and CLent*tion V Lto C "e > 'c^ n £^‘27“' 

bumpering information, and the incorporation of impact coeffici^ntc f 
matariais for .hich data is not av.if.bl.. Alio f£ opf£ia ! ng «r sion' T 

lltil /“ 1 f nprOVed t0 P rovide ”<>« choices for the variable being 
IdpnHf^rf 88 "k 1 - SS finish a11 as Pects of the choices that are currently 

allow wicks^ther ^han^ro!) 1 & ( f ’ 8 x fini * h the heat PiP e diameter variable to 
wiCKs otner than grooves to be used in the desien) Tt i e .1 „„ u j _ 

situations . mater 141 8 llbrary to offer additional choices ‘for design 

m . Ib if h ° Ped thiS ““"“a 1 is sufficiently clear to allow the new user to 
level^as h US6 ° f the HEPSPARC code. Similar information of different detail 

or to fh. JiSTitsS Zfrllfealfyl^f^dlon ^.^ 1 
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CODE LISTING 


A copy of the Fortran Source Code or ’ : Executable Image 

can Be Reguested By Contacting . 


Karl W. Baker 

National Aeronautics and Space Administration 
Lewis Research Center 
Solar Dynamics and Thermal Systems Branch 
21000 Brookpark Road 
Cleveland, Ohio 44135 
(216) 433-5278 



6.1 LIST OF SUBROUTINES 


The subroutines currently used in the HEPSPARC program are listed below 

serves^gb * 8 J nopsis of the section’s purpose and working. This list 

S reference to the user for determination of the intended 

function of the various routines used. The subroutines are listed in 
alphabetical order for ease of use of the section. 


Subroutine or Purpose 

function name 


...PAT 


XXXXXX 


The ... stands for a set of three letters (LITilithium 
NAK: sodium/potassium, HELrHelium, TOL:toluene, POTipotassium 
SOD; sodium, VAT:water, CES:cesium, MER:mercury, MET;methane, 
E:freon-ll, NHE:Heptane f DOW: Dowtherm A t and AMM:Ammonia) to 
represent subroutines that determine the compatibility of the 
working fluids with the various containment envelope materials for 
the heat pipes and transport ducts. 

Subroutines LITH ( lithium) , HEL(helium), TOUL(toluene) 

POT (potassium), SOD( sodium), CES(cesium), WATER(water) , 

HG (mercury), METH(methane) , FREON ( Freon-A) , NHEP(n-heptane) , 
DOWTHM( Dowtherm- a) , AMMON ( ammonia ) , STAIN( stainless steel) 

TITA( titanium), BERYL( beryllium) , ALUM( aluminum) , MAGN (manesium- 
lithium) , NIOB(niobium) , IRON (carbon- steel) , NICKEL(nickel) 

COPPER ( copper ) , TUNGS ( tungsten) , INCONN ( Inconel ) , 

MOLY (molybdenum), and TANT ( tantalum) . are all used to establish 
temperature dependent properties for the various materials needed 
in the radiator design. 


ACHANGE 


ARMOR 


This subroutine determines the change in the heat pipe vapor Mach 
number and other properties between the evaporator and condenser 
regions when a different wick structure is used in the two areas. 

This routine controls the code structure for determining the 
required armor thickness of the heat pipes and transport ducts. 
Various areas, fluxes, and thickness are determined based on 
geometrical and modeling flags passed to the routine. 


BLOCK DATA 
EFFICIENCY 


This block contains the information needed for the determination 
of the fin and surface efficiencies. 


DATA ThlS d ! tS blOCk that contains information required to determine 
FLUID some of the thermal properties of lithium, potassium, and sodium 

for use in the heat pipe heat transfer calculations. 

BRENT Tb f S is part of the actual minimization algorithm. This 

subroutine finds the actual minimum mass value of the radiator 
after MNBRAK has bracketed the value. This routine was taken from 
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CAPLIM 


CAPOUT 

CHECK 

COMPF 

DEBRIS 

DIAVAR 


DINTER 

DPLIQ 

DPVAP 

DPVC 

DPVE 

EFFS 

FACT 


Reference 17 and modified for use here. 

This routine determines the capillary limit of the heat pipe in 
question. 

This routine establishes key words for the compatibility 
determination of the fluids and materials. 

This routine contains the control structure to determine the 
compatibility of the working fluids and container materials. 

This computes fitted thermal functions for use with the FLUID 
BLOCK DATA information. 


This establishes a guess at the maximum thickness of armor 
required when a space debris flux is used for a radiator design. 

This subroutine establishes the parameters of the heat pipe during 
the optimization design option using the heat pipe diameter as the 
parameter. More than just the diameter needs to be varied for 
this type of analysis, and this algorithm accounts for changes in 
the wick and internal structure of the heat pipe as the outside 
diameter is changed. 


This function supports the EFFS routine. It interpolates values 
of the fin and surface data and also does some error checking. 


This subroutine calculates the pressure drop experienced by the 
fluid contained in the wick structure of the heat pipe. It is 
used in the determination of the limits of the heat pipes. 


This routine contains the control structure for determination of 
the vapor pressure drop in the heat pipe evaporator and condenser, 


This subroutine calculates the pressure drop experienced by the 
vapor in the condenser section of the heat pipe. This 8 " 8e 
the determination of heat transfer and fluid dynamics of the heat 

pipe. 

This subroutine calculates the pressure drop experienced by the 
vapor in the evaporator section of the heat pipe. This 8 “ 8 * t 
the determination of heat transfer and fluid dynamics of the heat 

pipe. 


This routine calculates the fin and surface efficiency o t e ea 
pipe and fin combination based on material and geometric 
specifications . 


This routine is not currently used by the program but has not been 
removed from the source code. 
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FINDT 

FLUX 

FRIC 

FUNC 

GROOVE 

HEADER 

HEAT 

HP 

HP1 

HPM 

HTC 


This routine determines the saturation temperature of the fluid 

in d thp P ri r «. baS ? d °? it8 pre88ure within the heat pipe. It is used 
£,tel ° n ° f th ' he “ Ch.r.ct.ri.tic. of £? 

This subroutine establishes the actual flu* r-ir, n 

•p.ce debris, present for th. *““ 

This subroutine calculates the friction factor required for the 
detonation of the pressure drop of the fluid in the transport 

This subroutine contains the control structure for calculating fh. 

of i th! Ubr0Utined f temineS the hydrodynamic and thermal properties 
calculations? "“ k ° > * t VU1 b ' “" d 111 «" transfer"'" 1 ** 

on or in It ° f th# tran8 P ort duct that has heat pipes 

lK 1S », rOUtine contains the control structure for determination of 
the heat transfer of the system. Calculations are Jo^“ 

femn”” 1 ! 16 the squired number of heat pipes in a system as well as 

analysis^ 68 ^ Changlng to different heat pipe fluids during the 

hoff 8 ^ br °i? tine calculate s the heat transfer capabilities of the 
heat pipe from the transport duct to the outside of the heat live 
This is compared to the heat transfer determined from QSHP and ? 
temperatures are modified as needed to balance these val!^. 

^ 8 a r ^r, eStabll8heS the choices of Possible heat pipe fluids 
for a given temperature. This routine is used multiple times in 

o£e ^ SCtiV ? lnpUt option when the design requires more than 

of t^smern 6 W ° rklng ^ t0 C °" r ^ range ° f ta “ pa -tures 

Infonnation concerning the physical details of each heat pipe are 

struc^uJe to call Ur tU 8 int ® ractive input * Thia also has control 
characterize ^ * U ^ data 

This subroutine determines the convective heat transfer 
coefficient of the fluid within the transport duct. 
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INPUT 

INRDWT 

JOINT 

LIMITS 

MAIN 

MASS 

MAT 

MAT1 

MAT 2 
MNBRAK 

OPTIM 

OUTPUT 


This subroutine crri.s out the majority of the interactive input 
of the program as veil as control structure for the remaining 
interactive input. Also, establishment of the saved data file 

“ done Jr. a. well a. some incidental calculations and flag 

setting. 

This routine contains the control structure used to readand write 
minified incut data and wick data. The input, or read 
option of this program is employed when non- interactive input is 

used. 

This subroutine calculates the heat transfer that occurs tram what 

or outlet regions of the transport duct for each panel. 

. t he various limits of the heat pipes 

an^determine” if they are being violated in the radiator design. 


The 

Also, 


This is the main control structure of the entire program 

St: 

masses are determined. 

Thi. subroutine control, the det.rmin.tion of the properties for 

ttt dt p”:. end tr.n.port duct fluid, end rn.teri.ls. 

This determines control of the selection of Web oic.s for the 
radiator desian. It is used for transport duct and heat p P 
uttt Jd trmtr m.t.ri.1. in the inter.ctiv. input option of the 

program. 

j&s-zss.'rx 

pipes in the interactive input procedure. 

This routine is part of the actual minimization Thi8 

SSoSe brackets the minimum v.l»e for Brent’s method of 
optimization (BRENT). This routine was taken from Reference 
and modified for use here. 

This routine contains the control structure for 
correct limits for the optimization program for as many time 
(i.e. several different heat pipe types) as is going to 
required in the program. 


This subroutine prints .11 the fln.1 output of the program. 


Some 
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PROPTH 

QSHP 

SCREEN 

SINTER 

SOLVE 

SOMETH 

SUBSTR 

THERTB2 

USER 

VAPRES 

WICK 

WIKTHR 

WIKWAT 


v. lu « ^^SlTSSr^SL" provide'th^final P ° W ' r 
specifications of the radiator design. P 1 

3S utr&K 1 *,' of the he,t 

resistance in that .action of th. haM^?® * tr, " ,fer 

heat pipe t with d a t fin 1 attached h to t it r ba8ed r on r0,n the ° ut8ide of a 

rs=S£Sap^«i , =5a- 

parameters, survival requirements, and redundancies An IMSL 

SJS.S.-2 *• — °'- 

Datennines therm, 1 rasistanca of tha actual heat pipe .apor. 

str * ss th “ is «- 
pipe fluid” 6 deteraines the condensate atate point for th. heat 

material. P library. IMterial th,t iS *’'* U * b1 ' *» *h. bullt^n 

Iaia«ted°bv^h. detemin r “* r ' sist *»<=» to heat transfer 
g erated by the vapor when it evaporates and condenses. 

lt\l°T lne ° b ! :ains the P h 7 sica l and hydraulic characteristics of 
the heat pipe wicks during the interactive input. araCteri8tics of 

of\\e 8U t^rand 0n H ta H inS H the COntro1 8tru cture for determination 
vick structure. hydrodynamiC characteristics of the heat pipe 

well I^other^f™*” 6 ? th * maSS ° f the wick ln the heat pipe as 
well as other information needed to determine the heat pipe mass. 


81 



VKRDWT 


is ."ffarsari.-sss \£ c. 


is used. 
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